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!>. Abstract: We develop a comprehensive description of three flavor neutrino oscillations 

/^ inside the Earth in terms of neutrino oscillograms in the whole range of nadir angles and 

■^ for energies above 0.1 GeV. The effects of the 1-2 mass splitting and mixing as well 

^~~l the interference of the 1-2 and 1-3 modes of oscillations are quantified. The 1-2 mass 

^T splitting and mixing lead to the appearance, apart from the resonance MSW peaks, of the 

QQ parametric resonance peak for core-crossing trajectories at Ei^ ~ 0.2 GeV. We show that 

^^ the interference effects, in particular CP violation, have a domain structure with borders 

^ determined by the solar and atmospheric magic lines and the lines of the interference phase 

S^ condition. The dependence of the oscillograms on the Dirac CP-violating phase is studied. 

We show that for sin^ 2^13 < 0.1 the strongest dependence of the oscillograms on 6 is in 
the 1-2 and 1-3 resonance regions. 

Keywords: neutrino oscillations, matter effects, leptonic CP violation. 
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1. Introduction 

Substantial future progress in neutrino physics will be related to the long baseline ex- 
periments as well as studies of the cosmic and atmospheric neutrinos. These studies are 
expected to fill some of the outstanding gaps in our knowledge of neutrino properties, such 
as the value of the leptonic mixing angle ^13, the type of the neutrino mass hierarchy, the 
octant of the mixing angle ^23 and the size of the Dirac-type leptonic CP-violation. They 
are also expected to improve the accuracy of the determination of the already known param- 
eters, such as the mass squared differences Attt-I;^ and jArTT-gj^l and the mixing parameters 

923- 



Q\2 and sin^ 26*2 



The key element of these experiments is that neutrinos propagate long distances inside 
the Earth before reaching detectors, and therefore a careful analysis of the Earth matter 
effects on neutrino oscillations is necessary (see, e.g., references in [1, 2] as well as Refs. [3, 
4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]). In a previous publication [1] we 
have studied these matter effects by making use of "neutrino oscillograms" of the Earth as 
the main tool. These are the contours of constant oscillation probabilities in the plane of 
neutrino energy and nadir angle. The plots of this type were introduced by P. Lipari in 1998 
(unpublished) and then appeared in several publications [22, 23, 24, 25, 1]. The oscillograms 
exhibit a very rich structure with local and global maxima and minima, including the 
MSW [26, 27] resonance maxima in the mantle and core of the Earth and the parametric 
enhancement [28, 29, 4, 5, 30, 31, 32, 22, 33] ridges for core-crossing neutrino trajectories. 
It was shown in [1] that all these features, including the local and global minima and 
maxima as well as saddle points can be understood in terms of various realizations of just 
two conditions: the generalized amplitude and phase conditions. We refer the reader to 
Ref. [1] for details. It has been shown that these oscillograms are very useful for gaining 
an insight into the physics of neutrino oscillations in the Earth and should help plan the 
future experiments as well as interpret their data. 

The analysis in [1] was performed in the limit of vanishing "solar" mass squared split- 
ting Ar?i2i- While this approximation is quantitatively well justified at relatively high 
neutrino energies {Ei^ > 3 GeV), it is less satisfactory at lower energies and also misses 
some important 3-flavor features of neutrino oscillations, most notably CP violation. In 
the present paper we extend the study of [1] to the case Am2i 7^ and consider 3-flavor 
effects in neutrino oscillations in the Earth, with the emphasis on the effects of the CP 
violating phase 6. We explore in detail 

• the effects of the 1-2 mixing and splitting on the oscillation probabilities, 

• the interference of the 1-2 and 1-3 (i.e., "solar" and "atmospheric") amplitudes, 

• the effects of and the sensitivity to the Dirac-type CP-violating phase 6, 

• the dependence of the oscillation probabilities on the neutrino mass hierarchy. 

We perform numerical calculations of oscillation probabilities and also develop a simple 
analytic approach to interpretation of the obtained results based on lines of three types 
in the neutrino energy - nadir angle plane. These are the solar and atmospheric "magic 
lines", i.e., the lines on which respectively the solar or atmospheric contributions to the 
transition amplitude approximately vanish, and the interference phase lines. Construction 
of these curves allows one to identify the regions in the experimental parameter space that 
are most sensitive to the effects of non-vanishing phase 6. Our results can be useful for 
planning experiments with atmospheric and accelerator neutrinos, as well as neutrinos of 
cosmic origin. 

In the present paper, as well as in [1] , we confine our consideration to the study of the 
oscillations probabilities. Accurate predictions for the event numbers and sensitivities of 
future experiments can only be done when an information on the corresponding detection 
efficiencies and systematic errors becomes available. Still, some general statements can be 



made and estimates done even in the absence of such an information; we plan to present 
the corresponding analysis in a future publication. 

Three-flavor effects in neutrino oscillations in the Earth have been considered in the 
past [34, 35, 36, 7, 23, 37, 38, 39, 40, 41, 11, 42, 43, 44, 45, 24, 46, 47, 48, 17, 49, 50, 51, 
52, 53, 54, 55]. The main new results of the present paper are the analysis of the neutrino 
oscillations in the Earth in terms of the aforementioned three sets of curves, and use of the 
neutrino oscillograms for a detailed study of the domains of the parameter space that are 
most sensitive to the value of CP-violating phase 6. 

The paper is organized as follows. In Sec. 2 we give general 3-flavor expressions for 
the oscillation probabilities in matter. We present exact analytic results for matter with 
constant density and introduce a factorization approximation. In Sec. 3 present the neu- 
trino oscillograms for different oscillation channels and discuss the effects of non-vanishing 
mixing and splitting in the 1-2 sector on these oscillograms. We consider features of the 
oscillograms for the inverted mass hierarchy. Sec. 4 contains the discussion of the effects 
of the CP-violating phase 6 and their analysis in terms of the three sets of special lines. 
In this section we also discuss the sensitivity of the oscillation probabilities to the phase 
6 and its dependence on neutrino energy and nadir angle (baseline length). Conclusions 
follow in Sec. 5. 

2. Three-flavor neutrino osciflations in matter 

2.1 Evolution matrix and probabilities for symmetric profile 

We consider mixing of the three flavor neutrinos, Vf = {ve^v^^Vr) . The mixing matrix 
U, defined through Vf = Uvjm where Um = {1^1,1^2,1^3)^ is the vector of neutrino mass 
eigenstates, can be parametrized as 

U=U23hUi^I-5Ui2. (2.1) 

Here the matrices Uij = Uij{9ij) describe rotations in the ij-planes by the angles 6ij^ and 
li = diag(l, 1, e* ), where 5 is the Dirac-type CP-violating phase. 
Evolution of the system in matter is described by the equation 

duf _ fUM'^U^ 



dt 



[--2Er^^' "'• '■'■'' 



where E,^ is the neutrino energy and M = diag(0, A?n2i, Ayrig^) is the diagonal matrix of 
neutrino mass squared differences with ^m^,i = tti'^ — rnf. V = diag(V^, 0, 0) is the matrix 
of matter-induced neutrino potentials with Ve = ^/2GpNe, Gp and A'^e being the Fermi 
constant and the electron number density, respectively. The evolution matrix S'(t, to) (the 
matrix of oscillation amplitudes) defined through v{t) = S{t, to) z^(to) satisfies the same 
Eq. (2.2) with the initial condition S{to,to) = 1. 

It is convenient to consider the evolution of the neutrino system in the propagation 
basis z> = (t'e, ^2) J^s) defined through the relation 

iyf = U23l5i> (2.3) 



with v = Ui3l-sUi2t^- According to Eqs. (2.2) and (2.1), the Hamiltonian H which 
describes the evolution of the neutrino vector of state P is 

H = ^t/i3 Uu M^ Ul Ul + V . (2.4) 

This Hamiltonian does not depend on the 2-3 mixing and the CP-violating phase. The 
dependence on these parameters appears when one projects the initial flavor state on the 
propagation basis and the final state back onto the original flavor basis. Explicitly, the 
Hamiltonian H reads 

A 2 Pl3 + •5l2 Cl3 ^A + 2Ve £'i./Am|i Sl2 Cl2 C13 Ta S13 Cl3(l - sfg ^a) 

H= ^^^^ ... cl^TA -si2Ci2Si3rA |. (2.5) 

\ ••• ••• cfg+sfs^fgrA 

Here ta = Am^i/^rri^-^, Cij = cosOij and Sij = sinOij. We introduce the evolution matrix 
(the matrix of transition and survival amplitudes) in the basis (z/g, 1^2, i^s) as 



^ee ^e2 ^e3 > 
S = I A^^ A22 ^23 
^3e ^32 ^33/ 



(2.6) 



This matrix satisfies the Schrodinger equation with the Hamiltonian H. Then, according 
to Eq. (2.3), the evolution matrix in the flavor basis S is 

S = USUI JJ = U23ls- (2.7) 

As follows immediately from the form of the Hamiltonian in Eq. (2.5), the amplitudes have 
the following hierarchy: 

^e2' ^2e ~ ^A, ^eg, % ~ Sl3, ^32> ^23 ~ '^IS'^A, (2.8) 

i.e., ^23 and A^2 are the smallest amplitudes. 

In terms of the propagation-basis amplitudes, according to Eqs. (2.6) and (2.7) (see 
also [43]), the matrix S is given by 



/ Aee C23^e2 + S23e"''^^e3 -S23^e2 + C23e"*^v4g3 

^3e ''23^22 ^ *2; 



C23^2e + ^23e*'^^3e 0^3^22 + 43^33 + ^MA' -523C23(^22 " ^Ss) + ^^r | ' 



\-S23^2e + C23e' ^3^ -S23C23(^22 " ^Ss) + ^^^ ■523^22 + ^23^33 + ^rr 

(2.9) 

where 

Kun = S23C23ie~'^ A23 + c'^ A-^^) , 



fifj. — ^23'-23Vc ^23 ^ "^ ^327 

(2.10) 



J^fiT — L23e ^23 *23'^ ^32 ' 



Krfj. = K^r{5 ^ -5, 2 ^ 3) , 

K — —K 
Notice that K^p {a, f3 = fi,T) are proportional to the small amplitudes A2^ and Ar^^- 



Since the Earth density profile is to a good approximation symmetric with respect to 
the midpoint of the neutrino trajectory and there is no fundamental CP- (and T-) violation 
in the propagation basis, the neutrino evolution is T-invariant in this basis, which yields [11] 

^2e = ^e2 ' ^3e = ^e3 ' ^32 = ^23- (2-11) 

Therefore for K^/s we obtain 

K^T = y423(cos2^23 cos(5 — ism6), 

K-T^ = yl23(cos2^23 cos 5 + isin(5), (2-12) 

K^^ = -Krr = ^23 ^^^ 26*23 COS 6 . 

Notice that the diagonal elements S^^ and Srr of the evolution matrix (2.9) depend on the 
CP phase only via cos 6, whereas 5ee does not depend on S at all. The latter is a consequence 
of our use of the standard parametrization (2.1) for the leptonic mixing matrix. 
The oscillation probabilities are expressed through the matrix elements of S as 

Pal3 = Pil^a^l'l3) = \Sf3a\'^ with a, /3 = e, /LX, T . (2.13) 

From Eqs. (2.13), (2.9), (2.11) and (2.12) one finds for the probabilities P^f^ 

P,^e = 43\AJ^ + sh\AJ^ + 2 S23 C23 Re(e-*'^^*2^,3) , (2.14) 

Pre = slslAJ^ + 4s\AJ^ - 2 S23 C23 Re{e~'^ Al^A^^) , (2.15) 

Pf^fi = |C23^22 + ■523^33+ 2 ■S23C23 COsSA^^l^ , (2.16) 

PfiT = |s23C23(^33 - ^22) + (cos 26*23 COS (5 + i sin (5) ^23^ • (2-17) 

The probabilities Pf3a are obtained from P^/? through the substitution 6 -^ —6:^ 

Pl3a = Paf3iS^-6). (2.18) 

All the results presented in this section are also valid for antineutrinos if one makes sub- 
stitutions 

S ^ -5, Aij ^ Aij, where Aij = Aij{V ^ -V). (2.19) 

Notice that the amplitudes of transitions (2.14) and (2.15) that involve fg are given by 
linear combinations of two propagation-basis amplitudes. The other amplitudes depend on 
three propagation-basis amplitudes. 

2.2 Eigenvalues 

We will refer to parameters 612 and Am2i as the 1-2 (or solar) sector, and to parameters ^13 
and ^m^i as the 1-3 (or atmospheric) sector. We use the same form of the mixing matrix in 



^Note that in a matter with an asymmetric density profile one would also have to substitute V ^ V , 
where V is the reverse profile corresponding to the interchanged positions of the neutrino source and 
detector [11]). This, in particular, means that one would have to distinguish Aij from Aji and Kpa from 



matter as in vacuum, with substitution 6ij -^ 9^. The eigenvalues of the Hamiltonian H^ 
are identified in such a way that H^ -^ Am'f-^^/2Ei^ when V ^ 0. For densities (energies) 
between the 1-2 and 1-3 MSW resonances we find for the normal mass hierarchy 

^ 2E^ ' ^ ' ^ 2E^ ^ ' 

For high densities (energies) which are far above the 1-3 resonance we have 

^m^Am|icf2 ^m^Am|ic|3 jjnr^y^ (2.21) 

^ 2E^ ' ^ 2E^, ' ^ ^ ' 

In a constant density medium the oscillation phases equal 

2(j)fi = AHji L with AHji = Hf - Hf. (2.22) 

There are two independent frequencies in the problem, AH21 and AH^2- 

For antineutrinos there are no level crossings, and with the increase of density (energy) 
the eigenvalues have the following asymptotic limits: 

^ ' ^ 2Ey ' ^ 2Eu 

We will discuss the level crossing scheme for the inverted hierarchy in Sec. 3.7. 

2.3 Amplitudes in medium of constant density 

Oscillations in a matter of constant (but trajectory dependent) density layers give a good 
approximation to the results of exact numerical calculations for neutrino oscillations in 
the Earth. As we have shown in [1] they reproduce rather accurately all the features of 
the oscillograms (at least for Amli = 0). Thus, we can use the exact analytic results for 
constant-density matter to clarify various features of the numerical results for the realistic 
Earth density profile. This also allows one to obtain a parametric dependence of the prob- 
abilities, in particular parametric smallness of certain contributions to the probabilities. 
In what follows we present the results for one layer of constant density. We mark the 
corresponding amplitudes and probabilities with the superscript "est" . 

The exact formula for the v^ -^ v^ transition amplitude in matter of constant density 
is 



Sf^ = 2ie'^"^ 



Ur,U;:{* smcpl^, - e-^™i[/ISC7-3*sin<Ar2j , (2-24) 

where C/™- and (f)"l are the elements of mixing matrix and the oscillation half-phases in 
matter. Using the expressions for U"^ and U"^ in terms of the mixing angles in the standard 
parametrization, we can rewrite Eq. (2.24) as 



where 



5-* = cos eiTsAlf + sin ^s^e"*'^'" A^f , (2.25) 



^cst ^ _ • gi^™ p^g ^m gj^ 26/^2 sin (/)^i , (2.26) 

^^|t ^ _ . gi0S sin 2^^ [sin ^^^ e'"'^^' + cos^ ^^ sin (/>^J . (2.27) 



Here c/)^ = 0^ + 0^. Since to a good approximation 9'^ f« 6*23 and S'^ ^ 6 [56, 42], 
the amplitudes j4*^|* and jd'^^* can be identified with A^2 ^^'^ ^es ™ -^^- (^•-'-4) ^^"^ (2.15). 
According to Eq. (2.20), between the two MSW resonances (p^ ~ Am'^2L/^^u- Above the 
1-2 resonance sin^ 6^ ~ 1, so that the second term in (2.27) is very smaU: 

cos^ «r. - ^ (^) = J tan^ 20, (-^ j , (2.28) 

which is suppressed as 1/V'^. Here E^2 ^ cos 2^12 A?7i2;^/2y is the 1-2 resonance energy. 
Consequently, up to the phase factor, the amplitude ^^g is reduced to the 2zv form with 
the parameters (^13,^32)- Notice that above the 1-3 resonance 0^ -^ /\m\2L/4:Ep. 
Similarly, for the v^ -^ v^ amplitude we obtain 

Sf^ = 1 + 2ie*^S|[/-|2sin<^™ _2ie"*^™2|C/;r3l'sin<A^2- (2-29) 

In terms of mixing angles, 

Tjm _ mm. „m m m iS"^ Tjm _ m m O 'if\\ 

and the amplitude can be rewritten as 

5;^^; = cos^ e™4f + sin^ 0™^|f + sin 2^^ cos (5"^^^ , (2.31) 

where 

4f = 1 + 2i e''^^^ sin^ 9-^2 sm </>"! , (2-32) 

A|f = 1 - 2i e-^-^sl cos^ 0™ sin 0™ + 2« e^"^™! sin^ ^™ cos^ 0™ sin (/.^^ , (2.33) 

A5f = i e^"^^! sin 9Tz sin 2^^ sin 0^i . (2.34) 

Notice that A'^^ has exactly the form of the corresponding 2u amplitude driven by the 
solar parameters. The amplitude A'?^ also coincides to a very good approximation with 
the corresponding 2u amplitude driven by the atmospheric parameters. The contribution 
of the solar mode to A'^^ is strongly suppressed by the factor sin^ ^13 in the region of 
1-2 resonance and by cos^ 9"^ above this resonance. Up to the overall factor sin 9"^!^ the 
amplitude A'^^ depends on the solar (1-2) phase only, and in general it contains double 
smallness: sin ^13 and the one related to the solar mode of oscillations. Again in the 
approximation 6*^ f» ^23 and 5'^ ^ 5 the amplitudes (2.32), (2.33) and (2.34) can be 
identified with the corresponding amplitudes in the propagation basis. 

For completeness, we present also the expression for the v^, -^ u^ amplitude: 

Sf^ = l + 2i e^"^?! cos^ 0™ cos2 0™ sin (/-^^ - 2i e'^'^S sin^ 0™ sin (/.^2- (2-35) 



It can be rewritten in the form convenient for use at low energies: 



5,"f = 1 + 2i e^-^^i cos^ 9T9 sin 6T, - 2i e'"^^^ sin^ e™ 



sin(^;5'l + cos"^l'^sin0^'l 



. (2.36) 



In formulas for the amplitudes one can interchange the phases in the exponents and 
sines using the following identity: 

sin (/.™e-^'^3i = sin ^^e"*'^"^ _ gin ,^^1, (2.37) 

where (p^^ = ^32 -|- 02i- Due to the level crossing phenomenon this phase interchange is 
convenient for studies of oscillation effects in different energy ranges. 



2.4 The factorization approximation 

The elements of the evolution matrix in the propagation basis S depend in general on 
Am2i, 612, Ar?7.|]^ and ^13. As follows immediately from the form of the Hamiltonian H in 
Eq. (2.5), in the limits Am^i -^ or/and S12 -^ the state 1)2 decouples from the rest of 
the system, and consequently, the amplitude A^2 vanishes. In this limit, A^^ (as well as 
ylgg and See) is reduced to a 2z^ amplitude which depends on the parameters Arn^^ and 
^13. We denote the latter by A a: 

^A(Amii, ^13) = A^3(Amii = 0) . (2.38) 

It is this amplitude that has been studied in our previous paper [1]; the corresponding 
probability equals Pa = I^Ap- 

In the limit S13 -^ the state P3 decouples and the amplitude A^t^ vanishes. At the 
same time, the amplitude A^2 reduces to a 2v amplitude depending on the parameters of 
the 1-2 sector, IS.m'^^ and 612. Denoting this amplitude by As we have 

As{Aml^,ei2) = A^2(% = 0) ■ (2.39) 

We will also use the notation Pg = |^5p. 

This consideration implies that to the leading non-trivial order in the small parameters 
si3 and TA the amplitudes A^^ ^'^^ ^2e below the 1-3 resonance depend only on the "solar" 
parameters, whereas the amplitudes A^^ and A^^ above the 1-2 resonance depend only on 
the "atmospheric" parameters, i.e.: 

Kl - Ae - As{Aml„eu) , E, < Eg , 

A^3 ~ ^3^ ~ AA{Aml, , ^13) , E, > Eg . 

In what follows we will call the approximate equalities in Eq. (2.40) the factorization 
approximation, since the dependence of the interference terms in the probabilities, Ai A^^, 
on solar and atmospheric parameters factorizes in this approximation. 

An additional insight into the factorization approximation can be obtained from the re- 
sults for constant-density matter obtained in the previous section. According to Eqs. (2.26) 
and (2.27), up to the phase factors one has 



where 



Aa -^ Af = sin 2^™ sin 0^, As -^ A^'f = sin 2^™ sin c^s , (2.41) 



ArrioiL 



(/)S = ^^\/{cos2ei2T2VE,/Aml,)^ +sm^2ei2, (2.42) 

<Pa = ^^'^ ^(cos 2^13 T 2VEjAml^Y + sin^ 2^13 ■ (2.43) 

Here the upper (lower) sign corresponds to neutrinos (antineutrinos). 

Due to the level crossing phenomenon the factorization approximation (2.40) is not 
valid in the {E^^Qj^) parameter space of the 1-3 resonance where 1-3 mixing is enhanced. 
Indeed, using Eqs. (2.20) and (2.21) we find that for the normal mass hierarchy 

0^2 ~0A for E.^Eg (2.44) 



and 



/,"! ^ J,_, f„^ I71 ^^ TPR 



4>s for E^^E^^. (2.45) 



However, the last formula is not correct in the energy region near the 1-3 resonance and 
above it due to the 1-3 level crossing. In particular, from (2.21) we obtain for Ey S> E^^ 

^^, « ^ (Ami,c?3 - ^mlA^) - ^^ - <\>\ , (2.46) 

where c/)^ is the 1-3 vacuum phase. At the same time, <ps ~ (p^- 

The results (2.26) and (2.27) allow us also to evaluate corrections to the factorization 
approximation. As follows from Eq. (2.26), beyond this approximation the amplitude A^2 
acquires an extra factor c^ and still depends on just one oscillation frequency, determined 
by the solar mass splitting Am2i- The factor cos 6^ decreases with increasing neutrino 
energy, approaching the value l/v2 at the 1-3 MSW resonance and further decreasing to 
cos ^55 ~ above this resonance. This leads to an additional suppression of the "solar" 
amplitude A^2 ™ matter at high energies, on top of the usual suppression due to the 
quenching of the mixing in the 1-2 sector. 

According to Eq. (2.27) the exact expression for A^^ differs from A^^ by an additional 
term depending on the 1-2 frequency. Above the 1-2 resonance we have cos^ 9^ <C 1 , and 
this additional term can be omitted: 



e3 



sin2^J^sin0^2e~''^^ 



(2.47) 



Hence, up to a phase factor the amplitude A^^ is also reduced to the standard two-neutrino 
form and depends on a single oscillation phase (j)^. We can rewrite the amplitudes (2.26) 
and (2.27) as 

A'^f = -i e^-^^i cos OTsAf {E^ < ^fg) , (2.48) 

g-i<^- ^cst ^ p^^2 ^m g-j^ 26*^ Af] , {E^2 ^Ey<. Ef^) . (2.49) 



e2 
e3 



In the general case of a matter of an arbitrary density profile, one can show, using 
simple power counting arguments, that the corrections to the factorization approximation 
for the amplitude A^2 ^^^ of order s^g, whereas the corrections to the "atmospheric" am- 
plitude ylg3 are of order ta [50] , in agreement with our consideration for constant density. 
The amplitude vl^g does not in general have a 2-flavor form, once the corrections to the 
factorization approximation are taken into account. 

3. Effects of 1-2 splitting and mixing 

3.1 Neutrino oscillograms of the Earth 

In this section we study the neutrino oscillograms of the Earth - contours of constant 
oscillation probabilities in the plane of neutrino nadir angle and energy - in the complete 2>v 
context (see Fig. 1). In our numerical calculations we take the matter density distribution 
inside the Earth as given by the PREM model [57]. Unless otherwise specified, we assume 



the normal neutrino mass hierarchy and use the following numerical values of neutrino 
parameters: Amj^ = 2.5 • lO'^ eV^, Aml^ = 8 • lO'^ eV^ tan 6*12 = 0.45 [2, 58, 59, 60]. 

Recall that the distance L that neutrinos propagate inside the Earth is related to the 
nadir angle of the neutrino trajectory Qi, through 

L = 2R^cos@,y, (3.1) 

where R^ = 6371 km is the Earth radius. The value @u = corresponds to vertically up- 
going neutrinos which travel along the Earth's diameter, whereas Q^ = ir/2 corresponds 
to horizontal neutrino trajectories. For < Q^ < 33° neutrinos cross both the mantle 
and the core of the Earth, whereas for larger values of the nadir angle they only cross the 
Earth's mantle. We call the corresponding parts of the oscillograms the core domain and 
the mantle domain, respectively. 

In the limit of vanishing Am^i the main features of the neutrino oscillations in the 
Earth are determined by the MSW resonance enhancement of neutrino oscillations and the 
parametric resonance enhancement in the core domain (see Ref. [1] and Fig. 1). The former 
leads to the appearance of the MSW resonance peaks: one in the mantle domain of the 
oscillogram at neutrino energies E^ ^ {Q — 7) GeV, and another one in the core domain at 
i?jy ~ (2 — 3) GeV. The latter (parametric enhancement) produces three parametric ridges 
in the core domain. In the notation of Ref. [1] they are ridges A, B, C. Ridge A spans 
the energy range from E'jy « 3 GeV at Q^ = to E^y ~ 6 GeV at the mantle-core border, 
where it merges with the mantle MSW resonance peak (see Fig. 1 middle and right panels). 
Ridge B starts from Ei, ^ 5 GeV at Qi, = and becomes nearly vertical for Q^, = 33°. The 
energy of the ridge C increases from E^ ^ 10 GeV to high energies when the nadir angle 
changes from Q^ = to Q^, = 33°. The location of the ridges weakly depends on the 1-3 
mixing angle. The ridges differ by the oscillation phase acquired in the core [1]. 

3.2 Oscillograms due to 1-2 mixing 

The oscillograms for ^13 = are presented in the left panels of Fig. 1. In this limit, 
according to Eqs. (2.9) and (2.39), 

l-Pee = \Aj'' = Ps. (3.2) 

Thus, shown are the contours of constant probability Ps- 

In the 2i/ case the oscillation probabilities depend on Am? and neutrino energy Ei, via 
the ratio Am?/E^. Therefore when oscillations are driven by the solar splitting, Amli, 
the oscillation pattern is shifted to smaller energies as compared to that due to Am'^^. 
Moreover, the 1-2 pattern differs from the pattern for the 1-3 mixing due to the large value 
of the 1-2 mixing. Indeed, as a consequence of the large 1-2 mixing the following new 
features appear. 

1. The oscillation length at the resonance is smaller than that for small mixing 

^™ = ^id^ ~ ^^ ' ^^-^^ 

where l^ is the vacuum oscillation length. 
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Figure 1: Neutrino oscillograms in the 3i/-mixing case. Shown are the contours of constant 
probability 1— Pee (upper panels) and 1— Pgg (lower panels) for Airi^i = 8x 10~^ eV^, tan^ 9i2 = 0.45 
and different values of ^13 . 
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2. The resonance energy is shifted to smaller values not only due to Am^i ^ ^^T-ii but 
also because of the factor cos 26*12 ~ 0.4: 

E^, = ^cos2e,,. (3.4) 

Here V is an effective (average) value of potential. 

3. The degree of adiabaticity is better than for the 1-3 mixing case, and therefore the 
oscillation probability in the mantle is determined by the potential near the surface of 
the Earth, V, averaged over the distance of the order of oscillation length. 

4. The oscillation length in matter, Im = 2it/AH2i, monotonically increases with energy, 
approaching in the limit E^ ^ oo the refraction length Iq = 2tx jV . Recall that for 
small mixings Im first increases with energy, reaches a maximum slightly above the 
resonance and then decreases. 

5. The jump of the mixing angle at the mantle-core boundary is small. Therefore, a 
sudden distortion of the oscillation patterns at Q^, = 33° is not as significant as it is 
for the small 1-3 mixing, especially below the 1-2 resonance energy. 

These features allow one to understand the structure of the oscillograms. In the mantle 
domain {Qy > 33°) the oscillation pattern for neutrinos is determined by the resonance 
enhancement of oscillations. There are three MSW resonance peaks above 0.1 GeV, which 
differ from each other by the value of the total oscillation phase. The outer peak (G,/ ~ 82°) 
corresponds to (/> ~ 7r/2, the deeper one at Q^ = 60°, to (^ ~ 37r/2, and the inner one 
{@u ~ 40°), to (f) = 5tt/2. Recall that such a large phase can be acquired due to a smaller 
resonance oscillation length (3.3) in comparison to the length in the 1-3 mixing case, in 
which only one peak with (p = tt/2 can be realized (see the upper parts of the panels in 
Fig. 1). The resonance energy is given by Eq. (3.4), and for the surface potential we find 

^f2« 0.12 GeV. (3.5) 

The ratio of the 1-2 and 1-3 resonance energies equals 



^12 _ f^^lA /^cos26'i2\ /^FlsA -. 1 



E^. \Aml-.J V cos 26*13/' \V12J 50 



(3.6) 



Here V13/V12 ~ 1.5 (see [1]), since for 1-3 oscillations we should take the average of potential 
along the whole trajectory. The estimate (3.5) is valid for two outer peaks. For the peak 
at Qu = 40°, V is larger, and accordingly, the resonance energy is slightly smaller. 

The width of the 1-2 resonance is larger and therefore the regions of sizable oscillation 
probability are more extended in the E^y direction as compared to the oscillations governed 
by the 1-3 mixing and splitting. 

The resonance energy in the core is Ef^ ~ 0.04 GeV. Therefore for Ei^ > (0.10 — 
0.15) GeV the 1-2 mixing in the core is substantially suppressed by matter. Furthermore, 
at the energies above the resonance energy in the mantle {E^ > 0.12 GeV) the mixing in the 
mantle is also suppressed. Therefore the peaks with Pmax ^ 1 at E^ > 0.12 GeV should be 
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Figure 2: Parametric resonance due to 1-2 mixing and splitting at the parametric peak at 0i, = 
25.1° and E^, = 0.2 GeV. Left panels: trajectory of the neutrino polarization vector in the flavor 
space, and projections of this trajectory onto three planes. Right panel: dependence of Ps (blue 
line) and of the oscillation phase (red line) on the distance along the neutrino trajectory. The stars 
along the blue line correspond to the beginning (1), the mantle-core boundary (2), the midpoint 
(3), the core-mantle boundary (4), and the end of the trajectory (5). 



due to the interplay of the core and mantle effects. In particular, the peak at Ei, ~ 0.2 GeV 
and Qi, ~ 25° is due to the parametric enhancement of the oscillations. It corresponds to 
the realization of the parametric resonance condition when the oscillation half-phases equal 
approximately <?!)mantic ~ 7r/2 and (/>core ~ 37r/2 (note that the total phase « 57r/2, and this 
parametric ridge is attached to the 57r/2 - MSW peak in the mantle domain). The spatial 
evolution of neutrinos in this peak and its graphical representation are shown in Fig. 2. 
In the left panel the blue lines present the trajectory of the neutrino polarization vector 
in the flavor space. Recall that in terms of the elements of the 2x2 evolution S matrix, 
5*11, 5*12, 5*21 = — 'S'*2 and 5*22 = S^i (where the last two equalities are the consequences of 
unitarity and hold in the basis where the effective Hamiltonian is traceless), the polarization 
vector is defined in the fiavor space as the vector with components sx = Re[S'J']^S'i2], 
sy = lm[SliSi2\, sz = |5'iip — 1/2. Then the Vg survival probability is given by Pee = 
|5ii|2 = SZ + 1/2 (see [1] for details). 

The ridge at Ei, ~ (0.12 — 0.15) GeV and Q^ ~ — 12° can also be considered as 
being due to the parametric resonance with a larger core phase: i;^core ~ 57r/2. However, 
for energies E^, < 0.15 GeV the in-matter mixing is nearly maximal both in the mantle 
and core, and the effect of the mantle-core density jump is small. As follows from Fig. 3 
(left panel), the positions of the MSW peaks and maxima of the parametric ridges are 
determined very well by the intersections of the lines that correspond to the coUinearity 
condition and the phase condition as in the case of oscillograms due to the 1-3 mixing [1] . 
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Figure 3: Interpretation of the oscillograms due to the 1-2 mixing and mass splitting, for neutrinos 
(left) and antineutrinos (right). Shown (in color) are contours of constant Ps as well as lines of 
various conditions that explain different structures of the oscillograms. The lines correspond to 
coUinearity condition (white), the generalized resonance condition (green), and the phase condition 
(black). 

Recall that in terms of the elements of the evolution matrices in the mantle, S"™, and in 
the core, S'^, these conditions read 



Re[5i"f5i^5i^] = (collinearity), Re[(5™' 5'=5")ii] = (phase) 



(3.7) 



In Fig. 3 shown are also the lines of the generalized resonance condition for symmetric 
density profile: Im[(S''^5™)ii] = (see [1] for details). The elements of the evolution 
matrix have been found by precise numerical computations. Apparently the ridges lie 
along the lines of the collinearity condition. 

There are also several intersections of the collinearity and the phase condition lines in 
the core domain of the antineutrino oscillogram (right panel) . This shows the existence of 
the parametric enhancement in the antineutrino channel. 

For El, > 0.3 GeV the oscillation length practically does not depend on the neutrino 
energy and is close to the refraction length, Iq = 27r/(v2G'i?A'^e). Therefore the lines of 
equal phases become nearly vertical. According to Fig. 3, the lines of zero oscillatory 
factor, (j) = irk with k = 1,2,3, that determine the so called "solar" magic lines (see [61]) 



14 



are at G^ ~ 54° in the mantle and at ~ 30° and 12° in the core domain. The oscillation 
probabilities become smaller than 0.5% above 4 GeV in the whole range of the nadir angles. 

At lower energies, E^ < 0.1 GeV, one finds a regular oscillatory pattern with ridges 
and valleys. The distortion of this pattern at the core-mantle boundary is rather weak due 
to the smallness of the difference between the mixing angles in the mantle and core. 

In the 2iJ context the oscillation probabilities depend on E^/Am'^ and the mixing 
angle. Therefore with the increase of the mixing angle the oscillatory pattern obtained for 
the 1-3 mixing (upper parts of the oscillograms) should continuously transform into the 
pattern due to the 1-2 mixing (apart from the trivial shift of energy). We find that with 
increasing 6 the parametric ridge A transforms first into the MSW peak in the mantle 
(1-3 mixing) and then to the outer MSW peak of the 1-2 pattern. The parametric ridge 
B transforms into the second MSW peak. The 1-3 core peak splits into two parts. One 
part transforms into the third MSW peak of the 1-2 pattern in the mantle. The second 
part merges with the parametric ridge C and appears as the parametric peak in the core 
domain at iiJjy = 0.2 GeV at large mixings. 

At high energies the patterns for neutrinos and antineutrinos are rather similar. 

Turning on the non-zero 1-3 mixing leads to the appearance of an 1-3 oscillation pattern 
at high energies and to the interference of the 1-2 and 1-3 oscillation modes. We will discuss 
two types of the interference. The first one is the interference of modes characterized by the 
solar and atmospheric frequencies. The corresponding interference terms in probabilities do 
not necessarily depend on the CP-violating phase. The second type yields the interference 
terms which depends on the CP-violating phase. We will call this the CP-interference. 

In the following subsections we will consider the effects of the inclusion of the 1-2 
mixing onto the 1-3 oscillatory pattern in different oscillation channels. We will compare 
the probabilities computed in the three-flavor (3j^) context and in the two-flavor (2z^) limit 
Am|^ -^ 0, 6*12 -^ 0. The 2z^ probabilities are computed as P(Am|^,sin^ 2^13). That is, 
we take a single mass splitting in the 2^ context to coincide with the largest mass splitting 
in the Sv case (the normal mass hierarchy). The oscillograms are computed for 6 = 0. 

3.3 Ug — V(, channel 

In Fig. 1 (upper panels) we show the probability 1— -Pee for three different values of sin^ 2^13: 
zero (left), small (middle) and relatively large (right). As follows from the figure, in the 
first approximation the oscillograms for non-zero values of ^13 appear as superposition of 
the 2u oscillation patterns produced by the 1-2 and 1-3 mixings with small interference 
effects. The smallness of the interference terms for this channel can be understood in the 
following way. According to Eq. (2.14), (2.15) and (2.18), the total probability of the Ue 
disappearance is 

1 — Pee = Pen + Per = Pe2 + -^eS- (3-8) 

When the 1-2 splitting is neglected, this probability reduces to Pa, studied in detail in [1]. 
Eq. (3.8) shows that for Am^i 7^ the z^e — ^^2 ^^"-^ ^e — i^^ transition probabilities simply add 
up in 1— Pee and no interference between the corresponding amplitudes (no CP-interference) 
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occurs. Correspondingly, the probability Pee does not depend on the CP-violating phase 
in the standard parametrization. It does not depend on the 2-3 mixing either. 

It follows from Eqs. (2.26) and (2.27) that P^g is driven only by the "solar" frequency 
(oc Am^i), at least in the constant-density approximation, whereas P^g depends both on the 
"atmospheric" parameters (Atji^^, ^13) and on the parameters of the 1-2 sector. Therefore 
the interference of the solar and atmospheric modes originates from P^g = l^ggp. The 
interference of 1-2 and 1-3 oscillation modes can be quantified (at least for the mantle 
region) using the expression (2.27) valid for constant density matter: 

1 - P^^f = cos^ 0^(1 + sin2 ff[l cot^ eT2)Ps'' + PT + cot 0r2 sin 20^ cos 0^^ A$f* Af . (3.9) 

Here A^^ and A^g^ are the 2z/ amplitudes defined in (2.41) (note that in terms of A'^^ 
and A^g^ this expression is valid in the energy range between the two resonances. In the 
other energy ranges one needs to take into account the level crossing phenomenon, which 
changes the labeling of the phases). The last term in (3.9) is due to the interference 
of the solar and atmospheric modes which comes from l^ggp. Apart from the product 
^A*^5*' ^^is term contains additional small factors. In the region of the 1-2 resonance, 
9^ ~ ^13, cot^j^ ~ 1 and therefore the interference term is suppressed by a small factor 
sin 2^13. In fact, all corrections to the main contribution, P5, are of the order of the 
small probability P^^ ~ sin^ ^13. In the region of the 1-3 resonance Pg^^ is small, and the 
interference term is suppressed by cot 9^ ~ r/^. Consequently, all the corrections to the 
dominant term P^^ are of the order of P^*^*. The interference term is further suppressed at 
the energies between the two resonances. Indeed, we can rewrite this term approximately 



as 



s'm29Y^sm29^cos(l)^iA'f^A'c^. As we will see in Sec. 4.5 for small ^13 the product 



2 0111 z,i/^2 °^^^ ^^^13 ^^° v^si^A 



sin 205^ sin 2^J^ has a minimum between the two resonances. Thus, the strong suppression 
of the effects of the 1-2 mixing in the fg — t'e channel is due to the absence of the CP- 
interference of the amplitudes A^2 and A^^. 

To further illustrate the effects of the 1-2 mixing and mass splitting at high energies we 
show in Fig. 4 the oscillograms for Pgg in the full three-flavor framework (left panels) and 
for the difference of probabilities with and without 1-2 mixing and splitting (right panels): 

APee = Pee " Pee = Pee " Peei^^l^ = 0). (3.10) 

Recall that we compute the two- flavor probability Pee taking IS.m? = Am|]^. 
In general, there are two contributions to APeg: 

APgg^APgi + APgi (3.11) 

which we will refer to as solar and atmospheric contributions. The solar contribution AP^ 
is proportional to the solar amplitude, it includes term P^^ and the interference of the 
amplitudes with the solar and atmospheric frequencies in P^g. The second contribution 
in (3.11) follows from a change of the atmospheric mode: the phase and the amplitude of 
oscillations due to non-zero 1-2 mass splitting and mixing: 

APg^ « AP^ = Pa-Pa- (3.12) 
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Figure 4: Oscillograms for the i^e — v^ channel. Shown are the contours of constant probability 
Pee (left) as well as constant difference APee of "iv and 2v probabilities (right), for neutrinos 
(upper panels) and antineutrinos (lower panels) . The oscillation parameters for iv probabilities are 
sin^ 26113 = 0.05, Am^^ = 8 x 10"^ eV^, tan^ 6*12 = 0.45 and 5 = Q. For the 2v probabilities we used 
Am^ = Arngj. 



Let us estimate these contributions using the results obtained for a matter of constant 
density. From Eq. (3.9) we have 



/^P^ « cos^ e'Tl Pt + cos eTo sin 2eX\ cos d)?, AT* A 



'13 ^u» VZI ^A ^S > 



(3.13) 



and, according to our consideration above, for high energies AP^ = 0{P^^) < sin^2^]^, 
which is below 0.002 in the 1-3 resonance region. For the atmospheric contribution we 
obtain 



AP^ ^ sin^ 2^^ sin2 <^™ - sin^ 29^3 sin^ 0™ 



(3.14) 
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Let us underline that due to the 1-2 level crossing, in the 3iy case the relevant atmospheric 
phase is cp'^ and not cp'^. In the lowest order 

sin2 2^i^ wsin2 2^J^(y ^Fi), (3.15) 

where 



V 
1 - sj^rA 
Then, using Eqs. (3.15) and (3.16), we obtain 



^1 - -. — :Tzr- (3-16) 



AP^ PS sin^ 2^i^ 



A(j) sin 2(/)^ - 2 ^^^^^ sin^ 0^ (cos 2i9i3 - x) xsfs^A 
where 



sin^ 2^5; 
sin^ 2(9i3 



(3.17) 



21/^^ 



x=— ^. (3.18) 

The first term in the brackets is proportional to the phase shift 

A<p = Ci - C2 - ^^ , (3.19) 

while the second one is due to the modification of the mixing angle. The doubly suppressed 
corrections, ~ Acpr^., are omitted. We find that for sin 2^13 = 0.05 both terms are of order 
0.02 — 0.03 in the region of the 1-3 resonance (in some regions of the nadir angles the second 
term can dominate). 

The difference AP can be rewritten as 



AP^ ^ 2Af sinie"^^ 



sin^ 29™" 
cos ^f^ Acj) - . „ /^ sin 0^2 (cos 2^13 - x) xsl^r^ 
sm 2c'i3 



(3.20) 



so that its proportionality to A^* becomes manifest. Consequently, AP vanishes along 
the atmospheric magic lines, A^^ = [61] (see Sec. 4.2 for details). 

Thus, for the zvg — t'e channel, AP^ dominates over APgg at high energies and therefore 
it describes the structure of the oscillograms for the probability differences. In particular, 
this explains the fact that oscillograms for APee repeat the structure of Pee with certain 
shift in energy. Partly the difference APge can be eliminated by modifying the 2v value 
of Am^Y-> taking Am? 7^ Arrig^. In certain energy range the phase shift effect can be 
eliminated. Notice also that AP^ is not proportional to A^ or even to the corresponding 
oscillatory factor. 

With the decrease of neutrino energy, the effect of the 1-2 mixing increases. Since 
AP oc E~^ and AP oc -E~^, at lower energies the effect of the solar contribution becomes 
important. The interference effects of different modes are suppressed in Pee, therefore the 
order of magnitude of the contributions to 1 — Pee due to nonzero 1-2 splitting can be 
readily estimated from the upper left panel of Fig. 1. 

The oscillograms for the antineutrino channel Ve -^ Ve are shown in the lower panels 
of Figs. 1 and 4. Now apart from the change of the sign of b one needs to take into account 
the change of the mixing and of the level crossing scheme. According to (2.23), the phases 
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become (/>a — (l^si ^^^^ 4'S — 'p2i- Therefore, the relevant oscillation phase in A'^^^ is (/)^, 
and the amplitude (2.27) can be rewritten for antineutrinos as 

^cst ^ -^14,^^ gjj^ 20^[sin (P'^^e-'^i^ - sin^ 0^ sin (/.™] , (3.21) 

where we used the phase exchange relation (2.37). At high energies sin^ ^^ ^ 1, and 
the second term in (3.21) is very small. For the normal mass hierarchy both 1-2 and 1-3 
mixings monotonically decrease with energy. The effect the of inclusion of the 1-2 mixing 
and splitting is illustrated in the lower panels of Fig. 4. Since for the antineutrino channel, 
in the absence of level crossing, the phase (f)^^ is relevant, the difference (p^^ — (j)^i is smaller 
than the corresponding difference in the neutrino channel. Thus, both IS.Pa and the total 
difference of the probabilities turn out to be smaller than they are in the neutrino channel. 
Furthermore, the solar contribution plays an important role now, determining the vertical 
structure of the oscillogram. 

3.4 Ue — '^fi and z^g — ^r channels 

The transition probability P^e = P^i^fj, -^ Ve) (see (2.14)) can be rewritten as 

P^e = 4|^,2l' + sisl^esl' + si^ 2023|^:2^e3l COs(./> " 5) , (3.22) 

where 

ct> = arg(.4*2^,3) . (3.23) 

Unlike 1 — Pgej this probability contains the interference term between A^2 ^"^^ ^e3 which 
depends on the CP-violation phase. Furthermore, this interference term is not suppressed 
by additional small factors as it happens in the t'e — T^e channel. From the unitarity of the 
matrix S in Eq. (2.6) we obtain for the product of amplitudes in the interference term 

^e2^e3 = "^32^33 ~ ^23^22 ' (3.24) 

i.e., A*'Ar^^ is proportional to the small amplitudes ^23 ^^"^ ^32- 

Since the amplitude A^2 is suppressed at high energies due to the smallness of the 1-2 
mixing in matter, in the lowest approximation we have 

P^e ~ Sin^ ^23|^e3P ~ sin^ ^231^^^. (3.25) 

The maximal value of the probability P^e — 523 — 0.5. 

In the constant density approximation (valid for the mantle domain) , using the ampli- 
tudes (2.27) we find at high energies 

P;f = slsPT + sin 2^23 cos 6TsAfAf cos(0ri + S) + ©(P^ )• (3-26) 

Comparing this expression with (3.22) we find that for energies between the two resonances 
the interference phase satisfies 

4> ^ -Ci • (3.27) 

The precise expression for the phase in the constant density approximation can be obtained 
from (p = arg(^^|**yl^|*) = arg(A^|). 
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Figure 5: Oscillograms for the z/^ — z/g channel. Shown are the contours of constant probability 
P^e (left) as well as constant difference AP^g of Sv and 2j/ probabilities (right), for neutrinos 
(upper panels) and antineutrinos (lower panels) . The oscillation parameters for Sv probabilities are 
sin^ 26*13 = 0.05, Amli = 8 x 10"^ eV^, tan^ 6*12 = 0.45 and 6 = 0. For the 2iy probabilities we used 
Am? = Am§i. In the right panels we also show the solar magic lines (black) and the lines which 
correspond to the condition (3.32). 

Fig. 5 illustrates the effect of inclusion of the 1-2 mixing and mass splitting on the 
'^fj. — T^e and v^ — i>e oscillation patterns for Ey >\ GeV. The left panels correspond to the 
2>v case, whereas the right ones show the oscillograms for the difference of probabilities 



AP^e = P^e - Pf^ei^ml^ = 0) 



(3.28) 



and, similarly, for Png and APng. As in the z/g— i^e case, this difference has two contributions: 



A P « A P"^ -I- A P 



(3.29) 
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The solar amplitude contribution AP£, is now dominated by the interference term and 
turns out to be much larger than in the Ve — ^e case. The atmospheric contribution is 
determined by the same AP = Pa — Pa as before, Eqs. (3.17) and (3.20), with the 
additional factor S23 ~ 0.5: 

AP^e ~ sl^AP^ + sin 2023 cos e^sAfAf cos(</.^i + 6) . (3.30) 

Now in the 1-3 resonance region, the interference is < 0.03 — 0.04, whereas the at- 
mospheric one < 0.015. Therefore, in the first approximation it is APfg that determines 
the structure of oscillograms for AP^e, whereas AP^ leads to some corrections to this 
structure. 

Using the factorization approximation we can write the solar contribution as 

AP;f, ^ cos923\As\ {2 sin 0231 Aa\ cos{(t>^^ + 6)+ cos 0231 As\) . (3.31) 



This formula corresponds to the one in Eq. (3.22) and is valid for the energies below the 
energy of the 1-3 resonance. In the right panel of Fig. 5 we show the lines of vanishing solar 
correction: AP^g = 0. This equality is satisfied where either \As\ = (the solar magic 
lines), or 

2 sin 023 1 ^A I cos (0?^i + 5) + cos 023 1 ^5 1 =0 (3.32) 



(white lines). If the second term is neglected the equality (3.32) splits into two conditions: 
Aa = (the atmospheric magic line, see Sec. 4.2) and cos(03]; + 5) = (the interference 
phase condition line). These lines form a grid which we will discuss in detail in Sec. 4. 
This grid describes rather well the domain structure of the oscillograms. Some deviations 
of the actual domain borders from the grid are related to the corrections from AP . 
Representing AP in (3.17) as AP = A'^^K, we can rewrite AP^e in the constant density 
approximation in the following form (up to corrections 0{Ps))'- 

AP^e ^ A^f [sl^K + sin 2023 cos 0^:^.4^"* cos((/.™ + 6)] . (3.33) 

Thus, with the atmospheric term taken into account, the lines A^^ = still determine zeros 
of the difference AP^g- However, the other lines of the condition AP^g = are shifted by 
the term ^23-?^ with respect to the solar magic and the interference phase lines. 

Notice that the corrections to the 2iy oscillograms are enhanced in the 1-3 resonance 
region, especially in the core domain at 3 — 4 GeV, where Aa is large due to the parametric 
resonance. Below 2 GeV the corrections increase with decreasing energy because AP^^ 



fie 



OC 



P^ OC 1/E^ and AP^g oc l/E^,. 

One qualitatively new feature of the oscillograms with Arn^i 7^ is that the absolute 
minima of P^e appear there as isolated points (see Fig. 5). In contrast, in the limit Am^i = 
the absolute minima of P^e (and P^g) never appear as isolated points in the oscillograms, 
but always form continuous curves (valleys of zero probability). This is unlike for the 
absolute maxima, such as the MSW mantle peak or the parametric resonance peak in 
the core region, where even in the limit of zero 1-2 splitting the value Pa = 1 (i.e., 
Pfie = S23PA = 523) is reached only at a few isolated points. This feature is a consequence 
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of the symmetry of the matter density profile of the Earth, and can be readily understood 
in the following way. 

The condition for the absolute minimum of the transition probability, P^g = 0, or 
Sefji = 0, where S is the evolution matrix, can be written as 

Re(5e^) = , Im(5e^) = , (3.34) 

and for a generic profile the absolute minima of P^g are found as the points where the curves 
corresponding to the two conditions in (3.34) intersect. However, due to the symmetry 
of the Earth's matter density profile, in the 2-flavor approximation (Am^i = 0) in the 
basis where the effective 2-flavor Hamiltonian of the system is traceless the transition 
amplitude is pure imaginary [11]. This means that the condition Re[S'e^] = is satisfied 
automatically for all values of E^y and Q^- Therefore, the zeros of P^e simply coincide with 
the contour curves Im[S'e^] = 0. For Am^-^ ^ this is no longer the case, because the 
2-flavor approximation does not apply. 

For antineutrinos (lower panels of Fig. 5) the corrections are again determined by the 
interference term with somewhat smaller atmospheric contribution. Therefore, one can 
see a domain structure with vertical lines. With the decrease of energy the maxima of 
the corrections inside the domains monotonically increase, since so do both the solar and 
atmospheric amplitudes. Notice also that for 5 = the positive corrections are larger than 
negative; the situation in the neutrino channel is opposite. 

As we pointed out in Sec. 2, for the inverse channel one has Pg^ = P^e{S -^ ~^)i where 
it has been taken into account that the Earth density profile is symmetric. 

According to Eqs. (2.14) and (2.15) the oscillation probabilities Pre and P^r can be 
obtained from the corresponding probabilities P^e and Pg^ through the substitution S23 -^ 
C23, C23 -^ —S23 [47]. The interference term has the opposite signs for channels including 
Ur as compared with those with f^, which can be obtained from the unitarity condition 
Pee + P^e + Pre = 1 and the fact that Pee does not depend on 5. 

3.5 Ufj, — I'fj, and I'r — t^t channels 

The v^ survival probability, P^^, for symmetric density profiles is given in Eq. (2.16). It 
can be rewritten as 

J^^l^l — l'-23^22 ^ *23^33l 

+ 2 sin 2^23 cos 5 Re [7113(023^22 + ^23^33)] + sin^ 26*23 cos^ 5\A^^\'^ . (3.35) 

Note that P^^ is an even function of 5. Since ^23 = C>{r/\si^) is a small quantity, one can 
to a very good approximation neglect the term ~ cos^ 5 in Eq. (2.16), which is proportional 
to |^23p- The term ~ cos^ 5 can only become important when the main terms in P^^ are 
small; however, this happens only in very small regions of the experimental parameter 
space. 

In the limit Amli ^ we have A22 = 1, ^23 ~ ^- Then, parametrizing the 33- 
amplitude as 

^-- ^ |.l33|e-^2^i = yr^T^e-^^^i (3.36) 
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we obtain from (3.35) 

P^^i^ml^ = 0) = 1 - sin2 2023 sin^ ,/>| - sI^Pa - sin^ 2^23 cos 2(/.|(l - ^l - Pa) . (3.37) 

If in addition ^13 = 0, then Pa = and P^^(Am2i = 0) is reduced to the standard 2u 
vacuum osciUation probabihty with (j)^- = 0^. 

It is easy to estimate the effect of the 1-2 mixing in the Umit ^13 = 0. In this case the 
eigenstate P3 decouples in the propagation basis, A^^ = 0, and the probability takes a very 
simple form 

Pnn = |c23^22 + ^23^33! • (3.38) 



Now 



and parametrizing A22 as 



A33 = e--^^3, 4 = ^> (3-39) 



A,, = |A,,|e-^^^r. = yr^T^e-^^^i , (3.40) 



we obtain 



P^^, = l- sin2 2023 sin2(,/)?- - ,/.?i) - 4P5 - \ sm' 2023 cos 2{<t)% - cb^g) (l - ^/l - Ps) . 

(3.41) 
At high energies, Ps <^ 1, and consequently Eq. (3.41) becomes 

P,, = P2- _ PscUcls + sis cos 2(4 - 0?!)] , (3.42) 

where P^f^ = 1 — sin^ 2023 sin^((/)?^ — 0^). In the constant density approximation and above 
the 1-2 resonance we have from (2.32) (j)^ = (p'!^, and therefore (p'^ — (p'^ = (p^. Thus, 
the effect of the 1-2 mixing is reduced to a shift of the oscillation phase and small additive 
correction of order Pg. Larger corrections are expected for non-zero 1-3 mixing due to the 
interference of the 1-2 and 1-3 modes. 

In Fig. 6 we show the oscillograms for P^^ and P^^ (left upper and lower panels) 
and the differences AP^^ = P^^ - P^^(Ami^ = 0) and AP^^j = Pp_p_ - Pp_p_{l\m\^ = 0) 
(right panels). One can see in this figure the regular oscillatory pattern dominated by the 
vacuum Vp, -^ v-r oscillations with certain distortion in the region of 1-3 the resonances: the 
MSW resonance in the mantle and core and the parametric ridges in the core domain. The 
dominant effect of the 1-2 mixing is related to the phase shift of this main mode, which 
increases with decreasing energy: at low energies the corrections become of order 1. AP^^ 
follows to a large extend the structure of P^^. No domain structure appears here. All these 
features can be seen from the formulas (3.35), (3.37). In particular, according to (3.41), 
for zero 1-3 mixing we obtain 

AP^^ = sin^ 2023 [sm'(</'^3 " C2) " ^i"' '^l^ + ^^^s) ■ (3.43) 

An additional insight can be gained using the constant density approximation. Let 
us consider the case of maximal 2-3 mixing for which the figures have been plotted. For 
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Figure 6: Oscillograms for the i/^ — z/^ channel. Shown are the contours of constant probability 
P^fj, (left) as well as constant difference AP^^ of 3iy and 2i/ probabiUties (right), for neutrinos 
(upper panels) and antineutrinos (lower panels) . The oscillation parameters for 3iy probabilities are 
sin^ 26li3 = 0.05, Am^^ = 8 x 10"^ eV^, tan^ 6*12 = 0.45 and 6^0. For the 2iy probabilities we used 
Am^ — Amli. 

maximal mixing and zero phase 5 the interference term in the f^ — u^ channel coincides, 
up to the sign, with the one in the u^ — Ue channel (see Sec. 4), which has been estimated 
in the previous subsection. Therefore, let us now consider the first term in (3.35). In the 
limit Am2i — > we obtain for this term a very simple expression: 



Pf^^iAml, = 0) = 1 - (1 - sin^ ^^3) sin' C2 ■ 



(3.44) 



Here the phase (j)^ should be calculated in the 2z/ context with ii^™ = 0. 

The survival probabilities Prr and Pff can be obtained from the corresponding prob- 
abilities P^^ and Pfip_ through the substitution S23 — > C23, C23 — > — S23 [47]. 
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3.6 u^ - u^ 



channel 



The probabihty of u^ -^ v-r oscihations for symmetric matter density profiles is given in 
(2.17). It can be rewritten as 

P^r = T sin^ 26*231^122 - ^33!^ + sin 2^23 cos 26*23 cos b Re [(A^g - ^13)^23] 

-sin26'23sin5Im [71*2^^3] + (1 - sin^ 26*23 cos^ (5) 1^23 p . (3.45) 

The oscihations are mainly driven by Am|^ and the large mixing angle ^23- The amplitude 
depends on 5 through the terms proportional to cos b and sin 5, and therefore P^t contains 
both CP- and T-even and odd terms. Due to unitarity, all the information on i/^ -^ v-^ 
oscillations is contained in the already discussed probabilities P^^ and Pp,e- P^it = ^ — P^ifi — 
P^e- Furthermore, one can show that the 5-dependent interference terms proportional to 
sin 6 and cos 6 satisfy the following relation 
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(3.46) 



(see the next section for details). 

Notice that for the maximal 2-3 mixing and 6 = 0, the probability takes a very simple 
form 
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^33! 



In the limit Am^i 



we obtain 



P^riAml = 0) = ^|1 - A33I2 = 1(2- Pa- 2^1 -Pa 



where 1?^^ 



4,- -.., 4v- -- -V- -^cos2</.| 
was defined in (3.36). For small Pa the probability becomes 

1 

" '-"--^32, 



P^ri^m^i = 0) = 1 - -Pa sin 



(3.47) 



(3.48) 



(3.49) 



f^ channel 



which in turn reduces to the standard 2z/ probability for ^13 = 0. 

The oscillograms for the t'^ — Vj- channel are very similar to those for the t'^ 
plotted in Fig. 6. To a large extent they are complementary in the sense that the corre- 
sponding minima and maxima are interchanged. The P^r oscillograms exhibit the vacuum 
oscillations pattern everywhere apart from the region Ey c^ 2> — 12 GeV. In this region the 
pattern is distorted by the 1-3 level crossing, as well as by the parametric enhancement of 
the oscillations in the 1-3 mode. This distortion is absent in the antineutrino channel. As 
in the Uf^ — u^ case, the difference of the probabilities, AP^,-, is dominated by the phase 
shift, and the corrections have the opposite sign compared to AP^^. 

Let us now present some results for constant density matter which will allow us to quan- 
tify the features described above. Using the expressions for the amplitudes in Eqs. (2.32) 
and (2.33), we find from Eq. (3.45) for E^ 3> Ey^ (the maximal 2-3 mixing and 5 = 0) 

2 
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e-^-^sl cos^ e'{^ sin 0^2 + [l " cos^ ^^(1 + sin^ 0^)] sin 0^1 



-i<t>Ti 



(3.50) 
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For energies above the 1-3 resonance one, cos^ 6^ -^ and Eq. (3.50) gives P^r ~ sin^ (/)^ 
with (/>2]; ~ Am'^iL/2E^. Therefore the corrections to P^r due to the vacuum 1-2 mixing 
and splitting are strongly suppressed. For energies between the two resonances, using the 
phase exchange relation (2.37), we obtain 

P;f « cos^ 0^ sin^ </," + - sin^ 20^^ sin (/)^i sin 0^i cos 0^2 + sin^ 6)^ sin^ 0™ . (3.51) 

For Am2i = we have 

P^fiAml = 0) « cos4 0™ sin2 </>- 

- - sin^ 26)^ sin (/.^^ sin 0^i cos ^2 + sm^ 0^ sin^ 0™ , (3.52) 

where the mixing angle 0^ and the phases should be calculated in the 2iy context. From 
the level crossing scheme for the normal mass hierarchy we obtain for energies below the 
1-3 resonance one that (f)^ ~ ~02i ^^^ 1^32 ~ ^31, and the difference between the 2z^ and 
Sv phases is proportional to Am\iL/2Ey. This difference increases with decreasing energy 
and at E'jy ~ 1 GeV it can be of order 7r/2. 

3.7 Inverted mass hierarchy 

Let us briefly comment on the features of the oscillograms for the inverted mass hierarchy, 
i.e., for Arrig-i^ < 0. The main change as compared to the normal hierarchy is due to the 
1-3 resonance structure which appears in the antineutrino channel now. Pulling out of the 
brackets in the Hamiltonian (2.5) the positive factor |Am|^|, we have to change the sign in 
front of all the terms without ta and V in the matrix (2.5). The general oscillation formulas 
we presented before are valid in this case, however the eigenvalues of the Hamiltonian and 
mixing angles in matter should be changed. The level crossing scheme is also modified. In 
the neutrino channel (where there is the 1-2 resonance only), in the limit of large energies 
we have 

^m^Arr4cf2 H^^v, H^ ^ ^^ . (3.53) 

In the antineutrino channel, below the 1-3 resonance the eigenvalues are 

Hr^V, H^^^^^, HT^^^, (3.54) 



2E,, 2E,, 



while above the 1-3 resonance, 



^m^Am|icf3 ^ Ami,cf2 H^^y. (3.55) 

^ 2E^, ' ^ 2E^ ' ^ ^ ' 

Here for illustration we will consider the oscillograms for the v^ — v^ channel only. The 
other channels can be analyzed in a similar way. In Fig. 7 we show the oscillograms for 
P^e (left panels), and for the difference of 2>v and 2v probabilities, AP^e (right panels). 

In the approximation of AroQ^i = the neutrino oscillograms for the inverted hierarchy 
coincide with the antineutrino oscillograms for the normal hierarchy, and vice-versa, pro- 
vided that Amg^ are taken to be the same in both cases. The inclusion of the 1-2 mixing 
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Figure 7: Oscillograms for the v^ — Vf^ channel in the case of the inverted mass hierarchy. Shown 
are the contours of constant probability P^e (left) as well as constant difference AP^g of 'iv and 2v 
probabilities (right), for neutrinos (upper panels) and antineutrinos (lower panels). The oscillation 
parameters for iv probabilities are sin^ 26*13 = 0.05, Am^i = 8 x 10^^ eV^, tan^ 012 — 0.45 and 
(5 = 0. For the 2iy probabilities we used Am^ — Am§]^. 



and mass splitting breaks this symmetry. However, at high energies where the corrections 
are smah, the correspondence "z^ inverted <-^ u normal" , "z^ inverted <-^ v normal" approx- 
imately holds. In our computations in the limit Am2i — > we have taken the remaining 
mass squared difference to be equal to the largest splitting in the Su context, Am^ = Am'^2- 
In the neutrino channel now there is only the 1-2 resonance and the height of the 
oscillation peaks of the probability increases with decreasing energy. The difference of Sv 
and 21/ probabilities is mainly due to the interference terms, as in the normal mass hierarchy 
case. This explains why one can clearly see a domain structure with vertical lines. The 
solar magic lines are shifted due to the contribution from AP . At high energies negative 
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corrections are larger in the absolute value than the positive ones (for 6 = 0). Notice also 
that the shape of domains for neutrinos in the case of inverted hierarchy is similar to that 
for antineutrinos in the normal mass hierarchy case. 

In the antineutrino channel the oscillation pattern is dominated by the MSW and 
parametric resonances. The difference of probabilities in the 3z/ and 21/ contexts is again 
due to the interference term oc AsAj^, with some corrections from APj^. According to the 
figure, the 1-2 mixing and splitting effect is enhanced in the regions of the MSW resonance 
peaks and along the parametric ridges, since A a is enhanced there. The oscillogram AP^e 
has a domain structure. The heights of the peaks are maximal in the resonance regions. 
In the mantle domain the heights reach minimum at i?,^ ~ 3 GeV and in the core domain 
at -E^ ~ 2 GeV, similarly to what we had for the CP-sensitivity peaks (Sec. 4). Due to 
the interplay of the interference terms and AP , for 6 = the corrections in the case of 
the inverted mass hierarchy are larger than those in the normal mass hierarchy case: e.g., 
in the region of the 1-3 resonance we have AP^e ~ ±(0.07 — 0.10). Below 3 — 4 GeV the 
positive corrections dominate. 

One can use the analytic formulas of Sec. 3 with appropriately changed phases and 
mixing angles to describe these results quantitatively. 

4. Effects of CP-violating phase 5 

In this section we consider in detail the properties of the CP-interference terms and in 
particular their dependence on the phase 6 in different channels. As can be seen from the 
expressions given in Sec. 2, the survival probability Pee does not depend on the CP-violating 
phase 6, both for oscillations in vacuum and in matter [34, 38]. This is a consequence of 
the facts that (i) 6 is rotated away by transforming to the propagation basis, and (ii) the 
probability Pee is not affected by this transformation. Note that for oscillations in vacuum 
or in matter with symmetric density profiles, the other two survival probabilities, P^^ and 
Prr, depend on 6 only through the terms proportional to cos 6 and cos 26 [43] since they are 
T-even quantities. In contrast to this, for oscillations in a matter whose density profile is 
not symmetric with respect to the midpoint of the neutrino trajectory, these probabilities 
acquire also terms proportional to sin 6 and sin 26. 

4.1 Interference and CP- violation 

The unitarity of the evolution matrix in the propagation basis (2.6) for symmetric density 
profiles gives 

"^23^52 + ^SS'^SS = ~^e2^e3 • (^-1) 

This relation allows one to write explicitly the 5-dependent terms for different oscillation 
channels as 

P^, = sin 2^23 {cos6Re[Al-^AJ+sm61m[Al^A£} , (4.2) 

P'^^ = sin 2023 cos6{- Re[Al^AJ - cos 2^23 M^U^m " ^22)] } > (4-3) 

P;!^ = sin 2023 {-sm61m[Al^A£ + cos6 cos2923Re[Al-^{A^^ - A^-^)]} . (4.4) 
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Here in Pr^ and Pr^T we have omitted small terms proportional to | ^23 1 ^ • The sum of these 
interference terms is zero. For maximal 2-3 mixing Eqs. (4.2), (4.3) and (4.4) reduce to 

Pi = cos 5 Re[Al^A£ + sin 6 lm[Al^AJ , (4.5) 

P'^^ = -cos6Re[Al^AJ, (4.6) 

P'^, = -sm5lm[Al^AJ. (4.7) 

The following consequences of these expressions can be useful for measurements of 6: if 
6 = 0, the probability Pnr does not contain the interference term and Pr^ = —P^e] if 



ference term and P^ = —P^e- 
Using the phase (p = ^^s{A*'A^n) defined in Eq. (3.23) we obtain in the general case 



6 = '7r/2, P^^ has no interference term and Pr^^ = —Pne- 



Pi= sin2^23Cos(</)-<5)|A,2^,3|, (4.8) 

P^^ = - sin 2(923 COS 5 cos (p \A^2A-^\ " -^23 , (4.9) 

P^^ = - sin 2^23 sin (5 sin </> |^e2^e3l + ^23 , (4.10) 

where 

D23 = ^ sin 4^23 cos (5 Re [^3(^33 - ^22)] (4-11) 

is proportional to the deviation of the 2-3 mixing from the maximal one. Notice that D23 
enters into P^ and P^ with opposite signs while P^ does not depend on D23 at all. D23 
is CP-even. It can be estimated using the constant density approximation as 

^23* ~ - ^ sin 4^23 cos 6 [cos (pf^c^P^Af A^'f - 2 sin ^^^ sin cp'^iAf] . (4.12) 

This expression shows that corrections to the term proportional to A'^^A^^^ in D23 are in 
general not small. 

For maximal 2-3 mixing one has 

Pi= cos(0- 5) 1^,2^,3!, (4.13) 

Pi = -cos5cos0|A^2^^3l, (4.14) 

Pi = -sin<5sin0|A^2^^3|. (4.15) 

The other (5-dependent terms in the probabilities, which are proportional to the square 

-|2 
23 1 ' 

Pfil = cos^ Jsin^ 2^231^23^ ' PfiT = (1 - cos^ (Jsin^ 2^23) 1^23^ • (4-16) 

Notice that the sum of these terms does not depend on 6 and equals |^23p- 

Let us present also the 5— dependent terms of the probabilities for other channels. As 
we have already mentioned, P^g = —Pne, and for the reverse channels, according to (2.18) 
we obtain 

Pf^ = sin 2^23 cos((/. + 6) \A^2A3\ ' (4-17) 

Pf^ = sin 2^23 sin 5 sine/. I ^g2^g3 1 + D23 ■ (4.18) 



of the small quantity j^g^P) ^-^^ 
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As was pointed out above, P-rr can be obtained from P^^ through the substitution S23 -^ 
C23) C23 -^ — S23: 

p4 = sin 2^23 cos 5 cos <j) \A^2K^\ - D23 ■ (4.19) 

For antineutrinos, according to (2.19), the probabiUties have the same form as the cor- 
responding probabiUties derived above with the changed sign of 5 and the amphtudes 
computed for the opposite sign of the potential. 

Thus, the 6 dependent terms in all channels are expressed in terms of two combinations 
of the propagation basis amplitudes, l^gg^ggj and D23. In the case of the maximal 2-3 
mixing (-D23 = 0), only the first combination enters the interference terms. Furthermore, 
for the channels involving electron neutrinos, ly^ — Ve and z^,- — z^e, only the first combination 
is relevant, even for the non-maximal 2-3 mixing. For generic values of 5 the CP-dependent 
terms in all the channels but v^, — v^ are of the same order. 

To assess the (5-dependent interference terms, one can consider the difference of the 
oscillation probabilities for two different values of the CP-phase: 

APCP(5) = p^^(5) _ p^^(<5,) . (4.20) 

In practice this would correspond to fit of the probability with the true value of the phase 
5 by the probability with some assumed value of the phase 5q. In Figs. 8, 10 and 11 we 
show some examples of the oscillograms for Ai-*]^ and AP^f for different values of 5 and 
(5o = 0°. Although the CP-oscillograms appear to have a complex structure, this structure 
can be readily understood in terms of the three grids of curves, which we consider next. 

4.2 "Magic" lines and interference phase lines 

Let us first consider the v^ -^ Ve oscillation probability, for which the equality 

APCP(5) ^ p^g(5) _ p^g(5o) = p5^(<5) _ ps^^i^s^) (4.21) 

is exact. The condition AP^e = is equivalent to 

l^eS^esI cos((/> -6) = \A^2Ki\ cos(0 - 5q) . (4.22) 

The same condition holds for the i^e — i^t channel. This equality is satisfied if at least one 
of the following three conditions is fulfilled 

(A) A^^iE,,e,)=0, 

(B) A^^{E,,e,)=0, (4.23) 

(C) (PiE,, e,) -6o = - [^{E,, e,) -6]+ 27rl . 

The last condition implies 

<j)iE^,e^) = i5 + 6o)/2 + TTl. (4.24) 
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Under conditions (A) and (B) the equality (4.22) is satisfied identically for all values of 6. 
In these cases the transition probability does not depend on CP-phase. It is also satisfied 
trivially if (5 = Jq + 27rn, since the true and assumed values of the phase coincide.^ 

In general, the conditions (A) and (B) can be satisfied at isolated points in the {@u, E^) 
plane only. Indeed, in order for A^^ to vanish, both its real and imaginary parts must be 
zero. The conditions Keying = and Im^^g = each define a set of curves in the 
{@u, Ey) plane, and the curves from one set can only intersect with those from the other at 
isolated points. The same applies to the condition A^2 = 0. In contrast to this, as we will 
discussed below, in the factorization approximation both the conditions A^2 — ^s = and 
ylg3 = A A = are fulfilled along certain curves in the oscillograms. This happens because 
the amplitudes As and Aa take a 2-flavor form. In the bases where the corresponding 2x2 
Hamiltonians are traceless, both A a and As are pure imaginary because of the symmetry 
of the Earth's density profile [11]. Therefore for As and Aa to be zero, it is enough to 
require that their imaginary parts vanish. So, instead of conditions (A) and (B) in (4.23) 
we will consider equalities As = and A a = 0. 

In the factorization approximation the conditions in (A), (B) and (C) define three sets 
of curves in the oscillograms (see Fig. 8), which play crucial role in understanding effects of 
CP violation. Along the lines determined by (A) and (B) the probabilities P^g) P^^ Pre and 
Per do not depend on the CP-phase. The other probabilities (as we will discuss later) only 
weakly depend on the phase along these lines. The lines shown in Fig. 8 were calculated 
in the factorization approximation, without assuming constant-density matter, by solving 
numerically the corresponding 2-flavor evolution problems with the PREM Earth density 
profile. 

In what follows we will consider these lines and their connection to the conditions in 
Eqs. (4.23) in turn. 

• Solar magic lines 

Let us discuss the condition As{Ej^,Qy) = 0. Notice that at A^2 ~ ^s = the "solar" 
contribution to the amplitudes of the v^ ^^ Ue and Ur ^^ t'e transitions vanishes, and 

P/ie ^ Pe/i ^ '^23l^e3l ' ''"'^ ^ '^''" ^ ''23l^e3l " V'^-^'-'J 

In Fig. 8 the condition As = determines nearly vertical lines at the values of the 
nadir angle @i, ~ 54°, 30° and 12°. This feature can be immediately understood using 
the constant density approximation. Indeed, according to (2.41) the condition As = is 
fulfilled when 

sm^s{Eu,@u) = 0- (4.26) 

As follows from (2.42), Eq. (4.26) is satisfied when 

L(e.)« , ^^'^^'^ n = l,2, ... (4.27) 

Am|i J {cos 26*12 =F 2VE^/Aml^)'^ + sin^ 26*12 



•^In Figs. 8, 10 and 11 we show the lines of the condition (C) which correspond to certain values of phases, 
which we denote by {S + (5o)flg. According to (4.13) the interference term vanishes if (f) — 5 + -k 12 + vrn. 
Comparing this last equality with (4.24), we find that along the lines of the condition (C) the interference 
term vanishes for value of phase 5 = —(5 + 5o)fig/2 + 7r/2 + -Kn. 
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Figure 8: Oscillograms for the difference of probabilities AP5e^((5) = Pfie{S) — Pf^ei^o) with Sq = 0°. 
Shown are the solar (black), atmospheric (white) and interference phase condition (cyan) curves. 
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Furthermore, at energies that are much higher than the solar MSW resonance energies in 
the mantle and in the core of the Earth, Ei^ > 0.5 GeV, the condition (4.27) becomes 

L{@u) ^ -^ . (4.28) 

Note that it is energy independent and determines the baselines for which the "solar" 
contribution to the probability vanishes. In the plane {Qu, Ey) it represents nearly vertical 
lines Qy ~ const. 

There are three solar magic lines which correspond to n = 1 (in the mantle domain) 
and n = 2,3 (in the core domain). The existence of a baseline (L ~ 7600 km) for which 
the probability of Ve ^^ '^^ oscillations in the Earth is approximately independent of the 
"solar" parameters (Am^-^^, 612) and of the CP-phase 6 was first pointed out in [62] and 
later discussed in a number of publications (see, e.g., Refs. [63, 64, 15, 65, 66, 61, 67]). This 
baseline was dubbed "magic" in [63] . The interpretation of this baseline as corresponding 
to vanishing "solar" amplitude A^2j according to Eq. (4.28) with n = 1, was given in [61]. 
In [61] it was also shown that for neutrino trajectories crossing the core of the Earth there 
exist two more solar "magic" baselines, corresponding to the oscillation phase equal vrn 
with n = 2 and 3, and the existence of the atmospheric "magic curves" was pointed out. 
The three solar "magic" baselines, for which the amplitude A^2 vanishes, can be clearly 
seen in the left panels of Fig. 1 . 

In the 1-3 resonance region and above, the factorization approximation becomes invalid 
since the angle 9^^ is large. As a result, A^2 ^^'^ ^s become substantially different. Indeed 
in the constant density approximation, Eq. (2.26), the equality A^2 = is satisfied when 

sin</>™(^^,e^) = 0. (4.29) 

For E^ <^ E^r^ we have (j)"!^ ~ (/)s. But for energies of the 1-3 resonance and above 
'P21 1^ 4's- III particular, for energies substantially above the resonance energy, cpl^ ~ 0^, 
and according to Eq. (2.46), 

4Ei,Tm 

Arrin 



He.) ^ -^^ ■ (4.30) 
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Thus, in the 1-3 resonance region, the condition sin (1)^ ^ transforms into sin (f)^ = 0, and 
the lines of condition A^2 ~ bend. This can be seen in Figs. 8, where the solar magic 
lines substantially deviate from the lines of AP5^ = 0. 

In the antineutrino channel no level crossing occurs and (p^ ~ (ps everywhere. 

• Atmospheric magic lines 

The atmospheric magic lines are determined by the condition Aa^E,^, Qi^) = 0. When the 
condition A^r^^E^, Q^) ~ Aa{Ei,, Qi,) = is satisfied, the "atmospheric" contribution to 
the amplitudes of u^ ^^ Uf, and Uj- ^^ Ug transitions vanishes. In this case, too, there are no 
effects of CP phase on the probabilities of oscillations involving i^g or z^e- 

The properties of atmospheric magic lines can be easily understood in the constant 
density approximation. As follows from Eq. (2.41), the condition Aa = is satisfied when 
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sin (pA = {(/)A = T^k, k = 1,2, . . .) or explicitly 



AE^^nk 



L(e,)« , ^^, k = l,2,... (4.31) 

Arnl^J{cos2eis^2VE^/Aml^y + sm^2ei3 



For energies which are not too close to the atmospheric MSW resonance energy, the con- 
dition (4.31) reduces to 

^''- \47rk±2VL{e,)\' ^^■'^^> 

which corresponds to the bent curves in the {@u, E^) plane. For very large energies, where 
Am^^/2E <^ V, the atmospheric lines approach the same vertical lines as the solar magic 
lines (4.28). 

Let us now consider the condition A^^{Ei^, Q^) = 0. In the constant density approxi- 
mation it gives, according to (2.27), 

sin (/.^2 = -e*"^™^ cos^ ^^ sin 0^i . (4.33) 

In turn, Eq. (4.33) implies two conditions which follow from the real and imaginary parts 
of the equality: 

sin ^^2 = - cos 0^1 cos^ 6*^ sin c/)^^ , (4.34) 

sin 0^1 sin (/.^i = 0. (4.35) 

Both conditions can be satisfied simultaneously only at certain points of the parameter 
space, which illustrates our general statement in the beginning of this subsection. However 
at high energies cos^ 0^ <^ 1, and the equality (4.33) approximately reduces to 

sin0^ = O. (4.36) 

Furthermore, since ^32 ~ <^a in the energy range above the 1-2 resonance, in this channel 
the factorization approximation works well. The atmospheric magic lines reproduce very 
well the lines of APJ^^ = 0. No interconnection of the atmospheric magic lines occurs. 

• The interference phase condition 

Let us now find the lines in the {@u, E^) plane which correspond to the interference phase 
condition (4.24). We shall call these lines the interference phase lines. Consider the condi- 
tion (4.24) in the factorization approximation. Although A^2 ^-^d A^^ are pure imaginary 
in the bases where their respective 2x2 Hamiltonians are traceless, their relative complex 
phase (p in any fixed basis is different from zero. It just equals the rotation phase between 
the aforementioned 2-flavor bases. In the constant density approximation (3.27) (p ~ —(PTi- 
In the energy range between the two resonances we have 

(4.37) 
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Figure 9: Grids of lines along which sin (pi = for different adiabatic phases (pi. Dashed lines 
correspond to (ps (black), cpA (red) and i/)^ (blue) obtained in 2j/ context. Solid lines correspond to 
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(4.39) 



AttI - 2{5 + 6o) 

This gives rather accurate description of the lines AP^e = below the 1-3 resonance. 

In the 3i/ framework (p"^ differs from 0^ in the 1-3 resonance region. Above the 1-3 
resonance cp'^i ^ (ps, and as we discussed before, the lines of (ps = const (Eq. (4.28)) depend 
on energy. Hence, the interference phase lines become nearly vertical with increasing energy, 
as can be seen in Fig. 9. 

Summarizing, the phase (p2i that enters into the amplitude A^2 nearly coincides with 
(ps below the 1-3 resonance. However, above the 1-3 resonance it approaches the vacuum 
phase (p^. In turn, the interference phase (p ~ cp^, which approximately coincides with (p^ 
between the 1-2 and 1-3 resonances, approaches (j)s with increase of energy above the 1-3 
resonance. The lines sin (ps = and sinKp^ = cross, whereas sincp^ = and sinijii^ = 
do not. Thus, in the region of the 1-3 resonance, where the factorization approximation is 
strongly broken in the 1-2 channel, interconnections of the lines occur: the contours of zero 
AP5^^ transform (interpolate) from the solar magic lines to the interference phase lines and 
vice versa. The interconnection is related to the level crossing phenomenon and reflects 
the level crossing scheme, it reflects the described change of the phases (p2i and (p^. To 
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illustrate this effect explicitly, we show in Fig. 9 the lines sin (ps = and sin (/)a = 0, which 
represent the magic lines (defined in the 2^ context), and the lines sin(/)2]; = 0, smcj)'^ = 
and sin(;^32 = which represent the conditions of vanishing amplitudes in the 3z/ context, 
that is, the contours of vanishing AP^^. Clearly, they do not coincide with the contours 
AP^^ = 0, since we have taken into account the phase factors only. This grid reproduces 
qualitatively well all the features of the lines shown in Fig. 8. 

There is no level crossing in the antineutrino channels for the normal mass hierarchy, 
and therefore there is no interconnection of the lines there. 

4.3 CP-phase domains for channels involving fg 

The solar and atmospheric "magic" lines and the interference phase curves allow a simple 
interpretation of the CP oscillograms. The solar (nearly vertical) and atmospheric (bent) 
curves divide the oscillograms into a set of domains, which are in turn divided by the 
grid of the interference phase curves into sub-domains (see Fig. 8). From these figures 
one can see that the interference phase curves are steeper than the atmospheric curves 
in the case of neutrinos and less steep than atmospheric curves for antineutrinos, in full 
agreement with Eqs. (4.32) and (4.39). This fact is related with the sign in Eq. (4.32). 
The probability difference AP^^{5) vanishes at the borders of these sub-domains: On the 
solar and atmospheric "magic" curves because the probabilities are 5-independent there, 
and on the interference phase curves because they correspond to cos{4> — 5) = cos((/> — 5q). 
The signs of the probability differences in the neighboring sub-domains are opposite, with 
the maxima of the difference in the central parts of sub-domains. 

As can be seen from Fig. 8, with changing the true (or assumed) values of 5, the solar 
and atmospheric grids remain unchanged, whereas the grid of the interference phase curves 
moves up or down, in accord with Eqs. (4.28), (4.32) and (4.39). The constant-density 
approximation results of Eqs. (4.28), (4.32) and (4.39) reproduce the main features of 
these curves quite well. 

As can be seen from the figures, the borders between the regions of the positive and 
negative CP-phase effect do not coincide exactly with the magic lines, especially in the 
regions of intersection of these lines. This indicates deviation from the factorization ap- 
proximation and is related to the level crossing phenomenon, as we have discussed in the 
previous subsection. 

4.4 CP-domains for channels not involving f^ 

As follows from our consideration in Sec. 4.1, the 5-dependent parts of the probabilities for 
the channels which do not contain Ug have more complicated structure than those which 
do. Apart from the term proportional the product of the amplitudes A^2-^*~ they contain 
contributions proportional to the deviation of the 2-3 mixing from the maximal one as well 
as terms P proportional to the square of the small amplitude ^23- I^ what follows we 
will neglect the latter. 

In the case of the maximal 2-3 mixing the (5-dependent parts of the probabilities are 
given in Eqs. (4.14) and (4.15). Furthermore, P^^ = —Pan and Pf^ = —P^t- From 
Eq. (4.14) it follows that for 623 = 45° the CP-oscillograms for the survival probability P^^ 
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Figure 10: Oscillograms for the difference of probabilities AP^^{S) ~ -P/j/j(<5) — Pfj.fj.{So) with 
Sq = 0° . Shown are the solar (black), atmospheric (white) and interference phase condition (cyan) 
curves. 
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can be interpreted in terms of the same grids of the solar and atmospheric "magic" curves 
As = and Aa = that we used for the analysis of the oscillograms for P^g and Pre (see 
Fig. 10). 

The key difference from the previous case is that now the CP-phase and the interfer- 
ence phase dependencies factorize. In the interference terms P/; and P^^ they appear as 
cos 5 cos (p. Therefore, the third grid describing oscillograms in these channels consists of 
the interference phase curves 

^^-^Ti = l + 7rn, (4.40) 

which, unlike the interference phase curves for the ly^ — Ue and Ur — t'e channels, do not 
depend on the values of 5 and 5o- The difference between the probabilities calculated with 
the true and assumed values of 5 (as well as the effects of 5 in general) vanishes on the curves 
of all three types and take maximum values in the central parts the domains delimited by 
these curves. The borders of the domains do not move with change of (5, and the only 
change that happens is that within each domain the probability varies proportionally to 
cos (5. 

The 5 dependent terms of the transition probabilities, Py^^ and P^„ are proportional 
to sin 5 sin 0. Therefore for these probabilities the interference phase condition reads 

,/.«0ri=™, (4.41) 

Again the borders of the domains do not depend on 5 and with changing 5 the interference 
terms vary as sin 5. 

Notice that the borders of domains are rather stable with respect to variations of the 
neutrino parameters. Since the grids are determined mainly by the phases, their depen- 
dence on the 1-3 mixing is weak. With the decrease of the 1-3 mixing angle in vacuum the 
grid lines becomes closer to the lines of vanishing AP^p. Since ^23 is experimentally known 
to be rather close to 45°, the discussed solar, atmospheric and interference phase curves 
give a rather good description of the CP-oscillograms for P^^ even when ^23 deviates from 
the maximal-mixing value (see Fig. 11).^ 

4.5 Sensitivity to the CP phase in the v^ — v^ channel 

Let us now identify the regions of the experimental parameters Qy and Ey for which the 
oscillation probabilities have maximal sensitivity to the CP phase 5. 

Consider the variation of the oscillation probabilities with varying 5 while all the other 
oscillation parameters are fixed. As follows from Eq. (4.8), the maximal variation of the 
probability P^g with 5 changing between 0° and 360° is 

AP^T" = max[P^e(5)] - min[P^e((5)] = 2 sin 2^23 I^e2^e3l • (4-42) 

(Note that a similar quantity was considered in [68]). This quantity is maximized when 
I^e2^e3l takes the maximum possible value. Let us now discuss the dependence of |^e2^e3l 



■^Actually, the correction due to 623 7^ 45° is of the order | cos 26123 J; 0.15. 
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Figure 11: Contour plots for the probability difference AP5jf (5) = P^f_i{5) — P^^{5o) with 5a = 0° 
and 5 — 180°. Upper panels: S23 — 0.4, lower panels: 533 — 0.6. Shown are the solar (black), 
atmospheric (white) and interference phase (cyan) curves. 

on the experimental parameters Qy and Ejj in the factorization approximation. For mantle- 
only crossing trajectories it is sufficient to use the constant-density matter factorization 
approximation, in which 



\A^^\ = sm.2B'^2 |sin05| , l^g^l = sin 26*^ | sin 0^1 • 

Thus, we have to find the maxima of the quantity 

A = sin 20]^ sin 20^^ | sin cj)^, sin < 



'13 



^^21 



^311 



(4.43) 



(4.44) 



with respect to Qy and Ey. To do this exactly is a rather complicated problem, and the 
result would be bulky and not easily tractable; fortunately, an approximate maximization 
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can be readily carried out by studying the energy dependence of the different factors in 
(4.44). 

First, recall that for energies Eu > 0.5 GeV the "solar" phase (ps is essentially energy 
independent. Therefore sin (ps can be considered a constant factor when maximizing A with 
respect to the neutrino energy. The maximum | sin (ps\ = ^ can be achieved by properly 
choosing the values of the baseline L. These are approximately equal to the central values 
of L (i.e., of cos 0jy) between the solar "magic" lines in the oscillograms of Fig. 8 (note 
that in the second band, due to the existence of the mantle-core boundary, the maximum 
is shifted from the center). 

Next, consider the extrema of the remaining factor, sin 20^ sm29^ \ sin(f)j^\, with re- 
spect to the neutrino energy for fixed L. To do this, we make use of the fact that the 
function 

f{E^) = sin 261^2 sin 20^ (4.45) 

varies with E^ significantly more slowly than sin 0^. If f{Ei,) were constant, the maxima of 
A with respect to Ei^ would coincide with the maxima of | sin 0yi | ; in reality, the (relatively) 
weak energy dependence of /(-Ej^) on Ei, leads to a slight shift of the exact maxima A from 
those of I sin (pA \ ■ The main effect of the energy dependence of f{Ei^) is actually to modulate 
the maxima of | sin(/)/i|. In other words, the maxima of AP™*^^ nearly coincide with the 
absolute maxima of | sin (ps sin (/)a\ (equal to 1), which are achieved by a proper choice of 
the values of the nadir angle and energy; the actual height of the local maxima of AP™*^^ 
is determined by the value of f{Ey) in these maxima. 

This is illustrated by Fig. 12. In the Q^ direction, the maxima of AP^^^ correspond 
to the maxima of | sinc/^^j. This fixes the values of the baseline L. The vertical "domain 
structure" {i.e., the structure in the energy direction) is due to the oscillatory dependence 
of I sin (pA I • For fixed L the energies En at which AP™^^ has maxima are approximately 
found from the condition | sin(/)^| = 1, or (pA ~ ^3iL = tt/2 + vrn. The heights of these 
local maxima are determined by f{En)- 

To find out which of the peaks of AP™^^ are the highest, one has to consider the 
maxima of the envelop function f{Ej^). We will do that here for the resonance channel 
(neutrinos); the analysis for antineutrinos can be easily performed along the same lines. 

The condition df /dEy = yields the third-order equation 

2x^ - 3a;2(cos 26*13 + ^a cos 26*12) 

+ x{l + 4rA cos 2^12 cos 26*13 + r\) - va cos 26*12 - r\ cos 26*13 = 0, (4.46) 

where x is defined in (3.18). For sin 2^13 < 1/9 it has three solutions, which correspond 
to two maxima of f{Ei,) and a minimum between the two maxima. Expanding in r^., we 
find for the position of the low-energy maximum 

xi ~ TA cos 26*12 + r\ sin^ 2^12 cos 26*13 , (4.47) 

or 

El ~ E^2 (1 + ^A sin 26*12 tan 26*12 cos 26*13) • (4.48) 
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Figure 12: Contour plots for the probability difference AP™f^ — maxP^e — minP^e for 5 varying 
between 0° and 360° and all the other oscillation parameters fixed. Also shown are the solar (black) 
and atmospheric (white) magic curves. We set Am2i = 8 x f 0^^ eV'^ and sin^ 26'i3 = 0.05. 

Here -Efl i^ ^^^ energy of the 1-2 resonance (which corresponds to x = r^ cos 2^12). Thus, 
the low-energy maximum of f{Ey) practically coincides with the low-energy MSW reso- 
nance, and the 1-3 mixing produces only a slight (©(r^)) upward shift of the position of 
the maximum. 

The minimum oi f{Ey) is given by 



X2 



while the second maximum is at 



COS26' 



13 



^A/l-9sin2 2%, 



X'i 



- cos 2^13 



^\^ 



9 sin^ 26^13 • 



^13 



(4.49) 



(4.50) 
cos 2^13. 



Recall that the MSW resonance energy in the 1-3 channel corresponds to x 

Consider the dependence of f{E^) on E^ for different values of the 1-3 mixing. For 
^13 — > we have X3 — > 1 and X2 -^ 1/2, i.e., E^ — > £^(3 and E2 -^ Ey^/2. Thus, in this 
limit the positions of the maxima of f{E^) coincide with the 1-2 and 1-3 MSW resonance 
energies, while the minimum is approximately in the middle between them. 

At maxima, the values of function / can be estimated as /(xi) ~ sin 2^13 and /(xi) ~ 
TA sin 26*12. 
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With increasing 1-3 mixing, the minimum of f{E^) shifts to larger energies: 

X2~^ + ^sin22^i3, (4.51) 

whereas the second maximum moves to lower energies: 

X3 PS 1 - - sin^ 26*13 ~ x^., cos^ 26*13 • (4.52) 



For sin^ 2^13 = 1/9 one has X2 = x^ = l/\/2, i.e., the minimum and the second maximum 
of f{Ey) merge. This just corresponds to the situation when the local minimum of the 
cubic function on the l.h.s. of Eq. (4.46) touches the x-axis. 

For sin 2^13 > 1/9 only the low-energy maximum of f{Ei,) persists. Thus, for these 
values of sin 2^13 the effect of the CP phase is maximal at the 1-2 resonance and decreases 
with increasing energy. For the other value, sin^ 2^13 = 0.05 (which in fact is not too 
small), we find X2 ~ 0.57 and X3 ~ 0.89. 

These results allow one to readily understand the oscillograms for AP™*^^ (Fig. 12), 
at least for mantle-only crossing neutrinos. For sin^ 2^13 = 0.05 the minimum of AP™^^ is 
situated at E'jy ~ 3 — 3.5 GeV {x ~ 0.57), which can be seen in the strip between the first 
and the second solar magic lines: the peaks increase in height both with energies decreasing 
and increasing from £^1/ ~ 3 — 3.5 GeV (in contrast to this, for sin 2^13 = 0.125 the height 
of the peaks would monotonically decrease with increasing neutrino energy) . Note that the 
situation is somewhat different in the strip between the Earth's surface {Qy = 90°) and the 
first solar "magic" line: the peak at -E;y ~ 5 GeV is actually lower than that at ~ 3 GeV 
because the baseline is relatively short, and at high energies the oscillation phase (pA is too 
small for the condition | sin(/>yi| = 1 to be satisfied. 

For core-crossing neutrino trajectories the const ant- density approximation is not in 
general applicable, and a different approach is necessary. The maxima of APf^^^ in that 
case can, in principle, be analyzed in the factorization approximation by making use of 
simple formulas for 2-flavor neutrino evolution in 3-layer matter density profiles obtained 
in [31]. In Fig. 12 one can see strong enhancement of the difference of amplitudes in the 
core domain at E^, ~ 2.5 GeV, which is apparently due to the mantle-core effect. 

4.6 Sensitivity to the CP phase in the u^ — v^ channel 

Let us now discuss the sensitivity of the survival probability P^^ to the phase 5. For other 
discussions of this issue see, e.5., [52, 54, 55]. From Eq. (2.16) one finds 

P^,^, = \C + Dz\'', (4.53) 

where 

C = C23A22 + ^23^33 > D = 2 S23 C23 ^23 ' ^ = COS 5 . (4.54) 

The maximum and minimum values of P^^ with varying 5 then correspond to the maximum 
and minimum of the modulus of the complex number C + D z when C and D are fixed 



42 



and z is allowed to vary between -1 and 1. A simple geometrical consideration then shows 
that for 

the minimum of P^^ corresponds to 5 = and maximum to (5 = 180° or vice-versa, so that 
for AP;^^ = max[P^^(5)] - min[P^^((5)] one finds 



(4.56) 



AP^7^ = \P^A^ = 180°) - P^^(<5 = 0°)| 

= 8 S23 C23 |Re[^!3(ci3A22 + ^^3^33)] I • 

If, on the contrary, 

2 S23 C23 1^23!' > |Re[43(4^22 + ^L^ss)] I ' (4-57) 

then for the maximal variation of P^^ with b one finds 

/|Re[^!=(ci3A22 + 5^3^30)] I V 

A-^MM ~ n \ h^S23C23|^23l ' (4.58) 

\ 1^23! / 

Note that yl23 is a small quantity, so that the condition (4.55) is satisfied in most of 
the parameter space. Exceptions are the regions where the "main" contribution to AP"?f^, 
i.e., \c22^22 + ■^23^33^' ^^ anomalously small, that is the regions along the magic lines. 
This is illustrated by Fig. 13, where we show the oscillograms for the maximum probability 
differences AP™f^ for 6 varying between 0° and 360°. The areas corresponding to the 
regions where the condition (4.55) is not satisfied occupy a rather small fraction of the 
parameter space. Moreover, APJ^^ is small in these regions, so that they correspond to 
low sensitivity to the effects of the CP phase. Note that these regions would never appear 
if one neglected the term of the order of |^23p ^^ ^^^ expression for P^^. 

For the regions where the condition (4.55) is satisfied, one can obtain a simple expres- 
sion for AP™'^^ in the limit ^23 = 45°. From Eqs. (4.56) and (4.1) one then finds 

AP-- = 2|Re(^,2^:3)|. (4-59) 

Thus, for ^23 = 45° the oscillograms for the maximum probability difference AP°?f^ are 
also governed by the solar and atmospheric "magic" curves, as well as by the interference 
phase curves (4.40), in full accord with our discussion in Sec. 4.2. 

5. Discussion and conclusions 

The main purpose of the present paper is to gain a physics insight into the complex pattern 
of full 3-flavor neutrino oscillations in the Earth. To this end, we presented a detailed de- 
scription of the three-flavor neutrino oscillation effects in the Earth in terms of the neutrino 
oscillograms, i.e., contours of equal oscillation probabilities or probability differences in the 
neutrino nadir angle - energy plane. 

We have found that for very small or vanishing 1-3 mixing the oscillation pattern ap- 
pears in the low energy region with large (maximal) transition probabilities below 0.3 GeV. 
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Figure 13: Contour plots for the probability difference AP™f^ — max P^^ — min P^^ for 6 varying 
between 0° and 360° and all the other oscillation parameters fixed. Also shown are the solar (black) 
and atmospheric (white) magic curves. We set Amli = 8 x 10^^ eV'^ and sin^ 26'i3 = 0.05. 

In the mantle domain the oscillation pattern consists of three MSW resonance peaks, which 
correspond to the oscillation phases 7r/2, 3it/2 and 5tt/2, and the parametric resonance 
ridge in the core domain at E^, ps 0.2 GeV and 0,^ ~ 28 — 30°. 

For non-zero 1-3 mixing the oscillograms consists of the low energy pattern, where the 
effect of the 1-2 mixing dominates, and the high energy pattern, determined mainly by the 
1-3 mixing and mass splitting, if the 1-3 mixing is not too small. The low energy pattern 
is modulated by the high frequency and small amplitude effect induced by the 1-3 mode, 
whereas the high energy structure is modulated by the low (refraction) frequency small 
amplitude effect due to the 1-2 mass splitting and mixing. 

We studied in detail the effect of the 1-2 mode on the oscillograms for energies Ei, > 
1 GeV. At these energies, if ^13 is not very small, the oscillation pattern is determined 
mainly by the 1-3 mixing and mass splitting, whereas the 1-2 mass splitting and mixing 
lead to small corrections. In the fg — ^e channel the interference of the 1-2 and 1-3 modes 
is strongly suppressed, and the effect of the 1-2 mixing is due to corrections to the 1-3 
mixing and the atmospheric phase. In the f^ — u^ channel the effect of the 1-2 mixing is 
dominated by the interference of the solar and atmospheric amplitudes. These corrections 
have a domain structure in the Ei, — O^, plane. In the f^ — f^ and f^ — Vr channels the effect 
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of the 1-2 mixing is essentially due to the corrections to the phase of the main (vacuum) 
oscillation mode. This consideration is important in discussions of the degeneracies of 
parameters in terms of the oscillograms. 

We studied the properties of the interference of the amplitudes A^2 ^^^ ^es ^^^i ^^ 
particular, the effects of CP-violation which are associated with this term. The structure 
of the interference term is simply and rather accurately determined in the factorization 
approximation, when one has A^r^ « AA{A.m'^i,6i3) and A^^ ^ As{Am2i,0i2)- This 
means that the dependence of the parameters of 1-2 sector and 1-3 sector factorizes in the 
interference term. This approximation does not work in the resonance regions. 

We showed that the interference term, and therefore CP-violation, exhibit a domain 
structure in the E^ — Q^ plane. The borders of the domains are determined by the grids 
of magic lines (solar and atmospheric) and by the lines of the interference phase condition. 
In the neighboring domains the sign of the CP effect is opposite. Beyond the factorization 
approximation the interconnections of the solar and phase condition Unes occur, which are 
related to the level crossing phenomenon. 

We studied the dependence of the probabilities on the CP-phase. The character of 
the dependence on the CP-phase is different for survival and transition probabilities. In 
the standard parametrization, the Vg survival probability does not depend on 6. For the 
survival channels t'^ — v^ and Ur — i^r as well as for the transition channel t'^ — I'r ah three 
sets of the lines - the borders of the domains - do not depend on 6. Within a given domain 
the interference term changes as oc cos 5 in the v^ — v^ channel, and as oc sin 6 in the t'^ — Ur 
channel. For the probabilities of transitions which involve Ue, i-e. Ve — t'/i and Ve — fr^ the 
interference phase condition depends on 5, so that with changing 6 the corresponding lines 
move. Therefore, the modification of the pattern of CP violation is determined by this 
motion of the interference phase lines. 

The phase 5 can affect significantly all the oscillation probabilities but u^ ^^ ^e- Thus, 
in principle, one can study the effects of leptonic CP violation by precision measurements of 
the energy and zenith angle dependence of these probabilities. We find that the strongest 
variation of the probability with 6 occurs in the region of the 1-2 resonance. For sin^ 2^13 < 
1/9 the local maximum appears in the 1-3 resonance region. The weakest variation is at 

Eu « O.S^fg. 

Many features of the oscillograms discussed in this paper are unobservable in the 
present and forthcoming experiments. The accelerator experiments cover only several pe- 
ripheral regions of the oscillograms, which correspond to large values of the nadir angles 
Qu > 77°. The large underwater and ice detectors have high energy thresholds E > 15 
GeV. Thus, the most interesting and structured regions of the oscillograms turn out to 
be uncovered. Measurements of the oscillograms with atmospheric neutrinos could be per- 
formed using multi-megaton water Cherenkov detectors, which will have sufficient statistics 
in the energy range E > 1 - 2 GeV, where the energy and angular resolutions are good 
enough [69, 70]. Detailed consideration of various features of the oscillograms performed 
in this paper can help develop methods which will enhance the sensitivity of future exper- 
iments to the CP-violating phase and other neutrino parameters. 
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